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Preface 


Nuclear thermal hydraulics is the application of thermofluid mechanics within the nuclear indus- 
try. Thermal hydraulic analysis is an important tool in addressing the global challenge to reduce 
the cost of advanced nuclear technologies. An improved predictive capability and understanding 
supports the development, optimisation and safety substantiation of nuclear power plants. 


This document is part of Nuclear Heat Transfer and Passive Cooling: Technical Volumes and Case 
Studies, a set of six technical volumes and four case studies providing information and guidance 
on aspects of nuclear thermal hydraulic analysis. This document set has been delivered by Frazer- 
Nash Consultancy, with support from a number of academic and industrial partners, as part of 
the UK Government Nuclear Innovation Programme: Advanced Reactor Design, funded by the 
Department for Business, Energy and Industrial Strategy (BEIS). 


Each technical volume outlines the technical challenges, latest analysis methods and future direc- 
tion for a specific area of nuclear thermal hydraulics. The case studies illustrate the use of a subset 
of these methods in representative nuclear industry examples. The document set is designed for 
technical users with some prior knowledge of thermofluid mechanics, who wish to know more about 
nuclear thermal hydraulics. 


The work promotes a consistent methodology for thermal hydraulic analysis of single-phase heat 
transfer and passive cooling, to inform the link between academic research and end-user needs, 
and to provide a high-quality, peer-reviewed document set suitable for use across the nuclear 
industry. 


The document set is not intended to be exhaustive or provide a set of standard engineering ‘guide- 
lines’ and it is strongly recommended that nuclear thermal hydraulic analyses are undertaken by 
Suitably Qualified and Experienced Personnel. 


The first edition of this document set has been authored by Frazer-Nash Consultancy, with the 
support of the individuals and organisations noted in each. Please acknowledge these documents 
in any work where they are used: 


Frazer-Nash Consultancy (2021) Nuclear Heat Transfer and Passive Cooling, 
Volume 6: Molten Salt Thermal Hydraulics. 
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1.1 


Introduction 


This technical volume describes the thermal hydraulic phenomena that are either unique to, or of 
particular relevance to Molten Salt Reactors (MSRs). It also provides an overview of the physics 
and chemistry of the salts typically used, and provides guidance on which modelling methods are 
best suited to make molten salt thermal hydraulic predictions, and how to use them. 


It is recommended that Volumes 2 (Convection, Radiation and Conjugate Heat Transfer) and 3 
(Natural Convection and Passive Cooling) are read before reading this volume, because it builds 
onand emphasises any areas of contrast with the descriptions of phenomena and methods applied 
in more conventional fluids. It is also relevant to consult Volume 4 (Confidence and Uncertainty) 
because of the comparatively large uncertainties that are found across a range of aspects in MSR 
analysis. 


The Use of Molten Salt in Nuclear Reactors 


The use of molten salts in a nuclear reactor has a number of advantages and features from a 
thermal hydraulic perspective: as a coolant, salts have high boiling points, so reactors do not need 
to be pressurised (unlike water cooled reactors); they can operate at high temperatures, giving 
high thermal efficiency and the possibility of cogeneration of hydrogen and process heat; they do 
not exhibit violent chemical reactivity with, for example, water (unlike sodium) and they have a high 
volumetric heat capacity, reducing plant size. 


However, it is not principally for thermal hydraulic reasons that salts receive attention as part of 
reactor design. Certain salts containing dissolved fissile materials can be used as a liquid fuel, for 
example in a fast or thermal spectrum breeder configuration, providing efficiency in fissile resource 
utilisation and waste minimisation’. Some of the broader nuclear aspects and technical challenges 
of MSRs? are described by LeBlanc (2010) and Serp et al. (2014). 


A wide variety of salt compositions have been proposed for use for a range of reactor purposes: 


* Primary coolant for designs with solid fuels. 


« Fissile and fertile material carriers for liquid fuelled designs, where salts of Uranium (U), 
Plutonium (Pu), Thorium (Th) or Minor Actinides (MAs) (for transmutation) are dissolved 
within the mixture. 


* Secondary and tertiary heat transport loops. 
Including by reusing existing spent fuel and by transmutation of long lived fission products, substantially reducing the volume 
and half-life of waste. 


In some references, the term MSR is restricted to refer to reactors where the fuel is dissolved in the salt; in this technical 
volume it is used in its broadest sense, when a molten salt is used as a working fluid in any role. 
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Introduction 


Molten salts are also used in ‘pyroprocessing’ of spent fuel and there are cases where a molten 
salt is used or proposed as a secondary power transport or thermal energy storage fluid for other 
high temperature reactor types, or in non-nuclear applications, such as Concentrating Solar Power 
(CSP). Usually two or three component salt mixtures are used to reduce the melting point to a 
practical temperature. 


Molten Salt Reactor Programmes 


The main experience of MSRs was gained at Oak Ridge National Laboratory (ORNL) under the 
Aircraft Reactor Experiment (ARE) and Molten Salt Reactor Experiment (MSRE) programmes (Fig- 
ure 1.1). 


These first MSRs were thermal spectrum liquid fuelled reactors (ARE was moderated with BeO; 
MSRE with graphite). The MSRE in particular demonstrated the feasibility of many aspects of 
the technology, operating from 1965 to 1969. It is described in contemporary and retrospective 
overviews by Haubenreich and Engel (1970), Rosenthal ef a/. (1972) and MacPherson (1985). 
Reporting for the MSRE programme was extensive and detailed, and this documentation is freely 
available from the United States Department of Energy (US DOE) Office of Scientific and Technical 
Information (OSTI)° and also from third-party websites‘. 


The US programme expected to build follow-on reactors from the MSRE (such as the Molten 
Salt Demonstration Reactor (MSDR), Bettis et a/., 1972 or Molten Salt Breeder Reactor (MSBR), 
Rosenthal et a/., 1972) and continued to create refined concept designs with improved safety and 
proliferation resistance (Gat, 1986, for example). There were also conceptual studies performed in 
a number of other countries, for example, in the UK at Winfrith (UKAEA, 1974), and at EIR (now 
PSI) in Switzerland (Taube, 1978, Taube and Heer, 1980). 


More recently, there has been a resurgence in interest, in some cases drawing substantially on 
these previous designs. MSRs are one of the six Gen IV reactor concepts®, and the subject of 
a number of development programmes (this list is not intended to be exhaustive — a selection of 
these reactors are characterised in more detail in Section 1.3): 


USA: The Advanced High Temperature Reactor (AHTR) and Fluoride-Salt-Cooled High-Temperature 


Reactor (FHR) (Figure 1.2) have been developed as concept designs for the last two decades 
as fluoride salt cooled, solid fuel reactors (Andreades et al/., 2016, Qualls et a/., 2017). 
More recently Kairos Power have begun to develop a variant (KP-FHR) towards commer- 
cial deployment. There are other active commercial developers, for example Flibe Energy 
are developing the Liquid-Fluoride Thorium Reactor (LFTR); TerraPower are developing the 
Molten Chloride Fast Reactor (MCFR); ThorCon International are developing the Thorium 
Converter (ThorCon) reactor and Elysium Industries are developing the Molten Chloride Salt 
Fast Breeder Reactor (MCSFR). 


Canada: Terrestrial Energy are currently developing the Integral Molten Salt Reactor (IMSR). 


UK: Moltex Energy are developing the Stable Salt Reactor —- Wasteburner (SSR—W). 


3 www.osti.gov 
4 moltensalt.org or energyfromthorium.com/pdf 
5 www.gen-4.org/gif 
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Figure 1.1: The MSRE layout and core (Rosenthal et al., 1972 and Rosenthal, 2009). 
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Europe: The Molten Salt Fast Reactor (MSFR) (Figure 1.3) is under development, most recently 
supported by the EU SAMOFAR and SAMOSAFERS* projects. Seaborg Technologies are 
developing the Compact Molten Salt Reactor (CMSR) and the Copenhagen Atomics Waste 
Burner is a small heavy water moderated design. 


China: Parallel development programmes are ongoing for Thorium Molten Salt Reactors (TMSRs) 
with solid and liquid fuels (Pioro, 2016, Chapter 14). There are plans for a small scale solid 
fuelled demonstrator to be developed soon (TMSR-SF‘1) with a larger power reactor to follow 
(TMSR-SF2), and a similar staged approach to a liquid fuelled variant (TMSR-LF1/2). 


Japan: Thorium Tech Solutions are developing the Thorium Molten-Salt Reactor (Th-MSR/FUJl). 


Russia: The MOlten Salt Actinide Recycler and Transforming (MOSART) reactor is being devel- 
oped by the Kurchatov Institute (Ignatiev et a/., 2014), and by collaboration with the SAMO- 
FAR project. 


High level (although not necessarily current) information on these reactors can be found in the IAEA 
Advanced Reactor Information System (ARIS) database’ or on the World Nuclear Association 
website®. 


In addition, a commonality in technology is starting to emerge with some designs of high power- 
density fusion reactor, which intend to use the same salt as some MSRs as a liquid cooling and 
tritium breeding blanket (Forsberg et al., 2019). 


Design Features of Molten Salt Cooled Reactors 


The candidate MSR designs exhibit a wide range of reactor sizes, shapes and flowpaths. The 
primary circuit flowpath is mostly a consequence of the choice of fuel and the neutron spectrum. 
The flow paths from the active core to the primary heat exchangers is via pipework in some designs. 
In others, the heat exchangers are immersed in the primary coolant (in a pool type reactor) or 
otherwise closely integrated (such as the MSFR, Figure 1.3). 


The choice of static or mobile fuel, neutron spectrum, salt type and emergency or Decay Heat 
Removal (DHR) approach distinguish the types of reactor under development — each is discussed 
below. The choice of these design features is closely interdependent, and will determine the reactor 
physics, thermal hydraulics, and material/structural performance assessments that are required for 
design and licensing. The features of some of the example reactors noted above (and the MSRE 
for comparison) are shown in Tables 1.1 and 1.2. 


Static or mobile fuel: Some reactor designs use solid fuels, mostly TRISO fuel, where encap- 
sulated fissile particles are embedded in graphite prismatic blocks, or in spheres forming a 
pebble bed. This is a very similar fuel design to that proposed for High Temperature Gas- 
cooled Reactors (HTGRs). The notable exception is the Moltex SSR—W design, which uses 
a conventional fuel assembly, similar to that used in Light Water Reactors (LWRs) or Liquid 
Metal-cooled Fast Reactors (LMFRs), except that the fuel is a molten salt within the (vented) 
tube of each fuel pin. 

® samosafer.eu 


” aris.iaea.org 
8 www.world-nuclear.org/information-library/current-and-future-generation/molten-salt-reactors.aspx 
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Figure 1.2: Layout of pebble bed FHR flow and heat transfer paths, showing the DRACS 
DHR system (Forsberg et al/., 2015). 
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Figure 1.3: Core of the MSFR (Cervi et a/., 2019). 
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Reactor Type Power (MWt) Shutdown/cooling 

MSRE Thermal, liquid fuel 8 Drain tank 

SSR-W Fast, static molten fuel 750 Reactor vessel cooling 

TMSR-SF1/2 = Thermal, pebble bed TRISO 10/384 Reactor vessel cooling/DRACS 

TMSR-LF1/2 Thermal breeder 2/373 Drain tank 

AHTR Thermal, prismatic TRISO 3400 DRACS 

KP-FHR Thermal, pebble bed TRISO 320 Reactor vessel cooling 

LFTR Thermal breeder, liquid fuel 600 Drain tank 

IMSR Thermal, liquid fuel 400 Reactor vessel cooling 

MSFR Fast breeder, liquid fuel 3000 Drain tank 


Table 1.1: Features of a selection of MSRs under development, and the MSRE for comparison. The 
information source for each reactor is given in the corresponding row of Table 1.2. 


Reactor Primary salt(s) Secondary salts(s) Reference 

MSRE "LiF —BeF.—ZrF,—UF, ’LiF—BeF» (66-34 mole %) Haubenreich and 
(65.0-29.1-5.0-0.9 mole %) Engel (1970) 

SSR-W (Na—U—Pu—MA)CIl, (fuel) NaNO3—KNO3 (‘Solar salt’, Moltex (2018) 
NaF—KF—ZrF, (coolant) tertiary) 

TMSR-SF LiF —BeF, (FLiBe, coolant) LiF -NaF—KF (FLiNak) Zhimin et al. (2016) 

TMSR-LF LiF —-BeF,—UF,—ThF, (fuel) NaF—BeF, Xu (2017) 

AHTR LiF —BeF, (FLiBe, coolant) KF—ZrF, (secondary and Lin et al. (2019) 

DRACS) 

KP-FHR LiF —BeF, (FLiBe, coolant) NaNO;—KNOs (intermediate) | Blandford et al. (2020) 

LFTR LiF,—BeF,—*"7UF,, (fuel) LiF —BeF, (FLiBe) EPRI (2015) 
LiF,—BeF,—**ThF, (blanket) 

IMSR UF, + fluorides Fluoride salt ARIS 

MSFR LiF—ThF,—5UF,—(Pu—MA)F, | LiF—NaF—KF (FLiNak) or Di Ronco et al. (2019) 


or LIF—ThF,—*“°UF, (fuel) 
LiF—ThF, (blanket) 


LiF—BeF> (FLiBe) 


Table 1.2: Salts used in a selection of MSRs under development, and the MSRE for comparison. 
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Other designs use molten/liquid fuels, where the fissile material is dissolved in the primary 
coolant and flows around the reactor primary circuit?. The heat source from fission is re- 
leased directly in the coolant (e.g. it is a heat-generating fluid). This has the advantage that 
no fuel assembly fabrication is necessary and fuel assemblies cannot ‘melt’ in an accident. 
In an emergency, the fuel salt can be drained into tanks of a sub-critical configuration by pas- 
sive equipment. The fuel salt can be continuously reprocessed while the reactor is operating, 
removing fission products and minimising excess reactivity in the core. More details on the 
fuel cycle and broader features of liquid fuelled MSRs can be found in Krepel et a/. (2014) 
and with an emphasis on the MSFR in Pioro (2016, Chapter 7). 


Neutron spectrum and fuel cycle: Both thermal and fast spectrum reactor designs exist, both of 
which can operate as breeders, including with a thorium fuel cycle or for the transmutation 
(‘burning’) of MAs from existing spent-fuel stockpiles. Thermal designs are typically moder- 
ated by graphite that is immersed in coolant or fuel salt. For a TRISO fuelled thermal reactor, 
the graphite is present in the fuel pebbles or prisms, as well as in reflectors at the edge of 
the core. 


Salt type: The choice of what salt composition to use depends on a range of considerations: 


¢ Melting point and vapour pressure. 


¢ Thermophysical properties (density, coefficient of thermal expansion, viscosity, thermal 
conductivity, heat capacity, infrared absorptivity). 


¢ Material compatibility, corrosion potential and susceptibility or resistance to contamina- 
tion by trace environmental species (such as oxygen, nitrogen or water). 


* Cost. 


¢ Solubility of fissile and fertile materials as well as fission products (including metals and 
gases). 


¢ Stability under irradiation and at high temperature. 
¢ Neutronic characteristics: moderating ability and absorption. 


* Neutron activation’? and tritium production"'. 


The coolant or carrier salts used in MSRs are generally chlorides and fluorides of alkali (Li, 
Na, K) or alkaline earth metals (Be). This position was arrived at during the initial reactor de- 
velopments described in Section 1.2. The precise composition of the salts used in a reactor 
design will, in some cases, be a proprietary and protected piece of information of a commer- 
cial technology developer. The composition is also not necessarily fixed during operation, 
because of the evolution of the presence of fission products, contamination by corrosion 
products, or material addition and removal through on-line reprocessing and fissile inventory 
control. 

8 The distinction of mobile vs. static fuel is made instead of liquid vs. solid because the ability of the fuel to move is the key 
feature, e.g. in the SSR—W, the fuel is liquid, but does not move except by circulation confined within its fuel pin. 

10 For example, in chloride salts °°CI (which has a natural abundance of 76%, with °”Cl comprising the remaining 24%) can be 
activated to °°C! a beta emitting isotope with a long half-life. This may require isotopic separation of °’Cl for fuels (Holcomb 
et al., 2011). 

11 For example, lithium based salts are a common choice, but the SLi + n reaction produces tritium, which in a high temperature 
reactor migrates through structural components, this can be minimised by enriching the salt to contain high purity of the ’Li 
isotope (>99.99% compared to the natural abundance of 92% 7Li and 8% SLi) — this increases the salt cost. The converse 


situation applies to lithium based salts used in fusion reactor blankets, which are trying to maximise tritium breeding, and 
so require °Li enriched salts (Forsberg et al., 2019). 
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Minimum operating temperatures are likely to be above the melting point of a salt mixture (by 
50 to 100 K) for example) to avoid solidification in the primary circuit. Therefore, the selec- 
tion of the composition of a salt mixture and the resulting operating temperature affects the 
selection of vessel materials according to their safe operating temperature, which is in-turn 
affected by corrosion, neutron embrittlement and creep. 


Emergency shutdown and cooling: During a postulated station blackout or loss of forced cool- 
ing flow event, there is a need to both shutdown the fission chain reaction and remove decay 
heat from the fission products in the reactor core. It is expected that if primary heat exchang- 
ers are operating when there is a loss of forced primary circuit flow, then natural circulation 
will be able to be established. However, where this is not possible, there are a range of 
options to assure shutdown and provide heat rejection to an ultimate heat-sink: 

Reactor vessel cooling: In common with other Gen IV reactors (such as several LMFR de- 
signs), some MSRs have an ‘always-on’ natural circulation flow path created by an 
annulus around the reactor vessel, either in an open-loop to atmosphere (such as an 
air-chimney), or in a closed loop via a heat exchanger to the environment". There is a 
constant heat-loss during normal operation, but when the vessel temperature increases 
in the event of an accident, the convective and radiative heat transfer are designed to 
increase to the point where it dissipates the decay heat load. In Sodium-cooled Fast 
Reactor (SFR) designs dating back to the 1980s, this type of system is referred to as a 
Reactor Vessel Auxiliary Cooling System (RVACS); similar designs applied to HTGRs 
are referred to as Reactor Cavity Cooling Systems (RCCSs). 

Drain tank: The fissile salt is drained under gravity into tanks with a subcritical geometry with 
an adequate reactivity margin — the draining can be initiated without operator interven- 
tion when elevated temperatures occur, via the melting of a ‘freeze valve’. The drain 
tanks may need to have a large thermal mass (e.g. be surrounded by water) and/or will 
be either actively or passively cooled, including by systems similar to those described 
for reactor vessels above. 

Direct Reactor Auxiliary Cooling System (DRACS): A heat exchanger is immersed in the pri- 
mary coolant in the reactor vessel, and connected by one or more natural circulation 
loops to an air-cooled heat exchanger in an air-chimney (Figure 1.2). 

The reactor vessel cooling and DRACS approaches still require shutdown to be ensured by 

neutron absorbing control rods or similar. Predicting the performance of these DHR systems 

with a suitably high-level of confidence for design and licensing will be one of the primary 
thermal hydraulic analysis tasks required. 


12 Neutron activation of Ar in the reactor vessel cooling air to 4TAr may make discharging it continuously an unacceptable 
radiological release. This is also a potential release route for tritium (Section 2.4.4). 
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2.1 


Technical Context 


Molten Salt Reactor Phenomena 


It is not uncommon for thermal hydraulics engineers to regard the prospect of performing analysis 
with molten salts with some trepidation; the unfamiliarity of them as materials creates an element of 
mystery around how they should be considered. Given the number of interlinked phenomena that 
an MSR exhibits, this is a reasonable initial view. However, for thermal hydraulic modelling aspects 
specifically, there is nothing particularly different about modelling molten salts compared to other 
working fluids, but how they are employed does require that any analysis is considered from the 
perspective of a fundamental understanding of the heat transfer processes. The main challenges 
and considerations in an analysis arise because of flows with low Reynolds number, the presence 
of participating thermal radiation (made more important by high operating temperatures), heat 
generation in the fluid, possible solidification, and the need to have salt-specific and temperature 
varying material properties, that are not necessarily known with high certainty. 


When considering the available guidance on the physical phenomena and their analysis for MSRs 
it is notable that there is significantly less detailed and consolidated documentation available from, 
for example, the IAEA and OECD NEA, compared to the scale of that available for LWRs and LM- 
FRs, as described in Volume 3 (Natural Convection and Passive Cooling) and Volume 5 (Liquid 
Metal Thermal Hydraulics). Other than the documentation from the MSRE programme, the most 
comprehensive general-purpose references are related to salt properties (Section 2.3). For me- 
chanical engineering analysis aspects, the state of the available literature reflects the diversity of 
configurations, and the lack of programmes that have progressed to detailed design, as well as the 
limited operating experience with demonstration reactors. This, along with more uncertainty and 
variability in fuel and coolant salt properties, leads to greater uncertainty in thermal hydraulic as- 
sessments compared to other advanced reactor types. This potentially requires additional efforts 
to establish the necessary confidence in any analysis performed, depending on the stage of design 
development and the margin available. 


The following reports and papers are a suggested primer giving a wide-ranging overview of the 
phenomena relevant to MSRs, their thermal hydraulic analysis, and the context into which that 
analysis will provide information: 


Williams and Clarno (2008): An ORNL review of salts and the reasons to choose between them. 
This is a Summary paper of detailed assessments of several different aspects of salt choice: 
thermochemical and thermophysical properties; neutronic performance (capture and moder- 
ation) and activation; material compatibility and corrosion, and material cost. 


Holcomb and Cetiner (2010): An ORNL overview of salt heat transport systems, with applica- 
tions not restricted to nuclear reactors, although it draws on the MSRE background knowl- 
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edge. This is a useful ‘signposting’ document to historical ORNL reports on aspects such as 
materials, pumps, valves, instrumentation, preheating, tritium handling and degassing, salt 
vapour pressure, chemistry and corrosion control, freezing, salt preparation and purification. 
The salt properties and heat transfer sections of this report have been largely superceded by 
more recent reports (Sections 2.2, 2.3 and 3.1). 


Lin et al. (2019): A group of experts from US National Laboratories, Universities and the US NRC 
produced a Phenomena Identification and Ranking Table (PIRT) for the AHTR, which can 
usefully serve as a basis for understanding the important aspects of solid fuelled MSRs. 


Diamond et al. (2018): Similarly to the above, a pre-PIRT was carried out for liquid fuelled MSRs', 
considering a range of potential design concepts, including fast and thermal spectrum vari- 
ants. 


Betzler et al. (2019): A description of the modelling and simulation functional needs, their interre- 
lation and prioritisation for the purposes of demonstrating the safety of MSRs for licensing. 


The last three references in particular provide descriptions of a wide range of phenomena that 
need to be considered in MSR analysis. The most relevant and important of these from the per- 
spective of thermal hydraulic modelling (as opposed to neutronics and reactor physics, materials, 
fuel and moderator performance, structural integrity or instrumentation and control) are described 
in the remainder of this section. Methods for modelling the relevant phenomena are discussed in 
Section 3. 


There are a number of features and phenomena that are related to the flow and heat transfer of 
molten salts: 


* Flowin the core, primary circuit, and DHR circuits could be laminar, or in the laminar-turbulent 
transition region in some designs, because of small flow paths and the moderately high 
viscosity of the salts. This could be in narrow fuel channels or gaps, pebble beds or heat 
exchangers. This contrasts with flows in most reactor designs with non-salt coolants which 
are usually fully turbulent for processes of interest. Modelling this will require pressure drop 
and heat transfer correlations that are specific to this flow regime. 


Uncertainty and temperature dependence in the value of the viscosity of salts will lead to 
uncertainty in the primary circuit pressure drop and in natural circulation flow rates, both of 
which are important characteristics. 


Compared to LWRs, low pressure operation means that energetic vessel or pipe break flows, 
depressurisation, and the associated two-phase flow phenomena are not possible, simplify- 
ing some aspects of the thermal hydraulic analysis of an MSR. 


¢ There are expected to be geometrical features (such as open plena) and reactor conditions 
where asymmetric or 3D flows occur. 


« Reactor designs with DRACS DHR systems include a ‘fluid diode’ in the core flowpath that 
includes the DRACS heat exchanger, which has a high flow resistance under normal oper- 
ation to limit thermal losses, but a low resistance when the flow direction is reversed under 
loss-of-flow conditions. 


1 Considered ‘pre-’ because the state of knowledge of candidate MSR designs and the selection of licensing basis events 
was not considered to be sufficiently mature. 
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* Thermal spectrum MSRs contain graphite, either as part of the fuel or as a moderator. The 
dimensions and properties (including heat capacity, thermal conductivity and coefficient of 
thermal expansion) of graphite change under irradiation, as is well studied in Advanced Gas- 
cooled Reactors (AGRs), and recently measured under high temperature conditions by Hei- 
jna et al. (2017). Changes in fuel prism or moderator channel geometry with core life will 
affect heat transfer, pressure drop and flow distribution in the core. Changes in graphite tem- 
perature will also affect the reactivity of a liquid fuelled core, because, by expanding, it can 
reduce the quantity of fissile material in the core. 


Aspects relevant to understanding and modelling these phenomena are described in Section 2.2. 
Further features and phenomena that relate to MSR design and operation are: 


* Temperature effects on reactivity are particularly important in MSRs, because there is a 
strong negative reactivity feedback from Doppler broadening with increasing fuel tempera- 
ture. For liquid fuelled reactors this coupling is stronger, because an increase in temperature 
also reduces the fluid density, including reducing the quantity of fissile material in the core 
in some designs. The effect of any reactivity increase changes the primary fluid temperature 
directly (it is a heat generating fluid), so the feedback is immediate. 


* Temperature feedback can be coupled to laminar-turbulent transition flows — where a transi- 
tion occurs, the salt mean temperature (and the shape of its profile in a channel) will change, 
and potentially induce a reactivity change. This change in reactivity may lead to a change in 
buoyancy forces, that could interact with the generation or suppression of turbulence, and so 
this may be a coupled effect. 


¢ For liquid fuelled reactors, the transport of Delayed Neutron Precursors (DNPs) is a unique 
feature, because the source of an important proportion of the neutrons required for criticality 
and reactivity control moves through the reactor in the interval between absorbing a neutron 
and fissioning. This alters the reactivity and kinetics of the core by affecting keg (effective 
neutron multiplication factor) and Ber (effective delayed neutron fraction): 


— The movement changes the location of reactivity in the core. 


— DNPs that flow out of the core emit their neutrons elsewhere in the primary circuit 
(exposing pumps, piping and heat exchangers to a neutron flux). 


— DNPs that travel around the whole primary circuit before decaying will then flow to the 
core entrance, creating a reactivity source. The complexity of this increases during tran- 
sients, such as a reactor start-up. This will require particular attention when assessing 
the dynamic stability of the reactor. 


¢ Similarly, for liquid fuelled reactors it is necessary to quantify the salt composition/inventory 
and mass transport throughout the primary circuit through time to account for fuel deple- 
tion (requiring the neutron flux distribution) and material addition, removal and deposition. 
This will be used to evaluate reactivity as well as corrosion (redox) potential, thermophysical 
properties of the salt and radiological source terms. 


¢ The level of uncertainty in aspects of MSR material properties and operation, and the in- 
terdependency of the reactor phenomena, mean that coupled multiphysics methods and 
Sensitivity Analysis (SA) and Uncertainty Quantification (UQ) will be necessary, particularly 
for neutronic feedback, chemistry and species tracking. 
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« Freezing and melting of coolant or fuel salts may play a role in normal operation (via ‘freeze 
valves’) or in fault conditions. ‘Freeze plugs’ are intended to melt at temperatures that occur 
during accidents, draining the salt from the core. Conversely, overcooling faults could freeze 
salt in pipework and heat exchangers, preventing flow. This is also a concern where salt 
vapour may condense on surfaces in cover gas systems and pools. 


« Dissolved gases can be present in molten salts, and their transport and removal needs to be 
understood. They can be generated in the form of tritium or fission products produced by the 
reactor, or they may be inert gases injected to remove these species, that are themselves 
entrained in the salt. Transport of soluble gases can create localised reactivity feedbacks due 
to the evolution of bubbles or high neutron absorption cross-section species such as ‘Xe. 


The aspects of these phenomena that pertain to developing the thermophysical properties nec- 
essary for thermal hydraulic modelling and analysis are discussed in Section 2.3. Additional phe- 
nomena relating to salt chemistry, nuclear effects and multiphase behaviour are described in Sec- 
tion 2.4. It is important to be aware of these, but, because they are not directly part of single-phase 
thermal hydraulics analysis, they are not covered further in Section 3. 


Many molten salts are not fully opaque or fully transparent to thermal radiation (unlike, for example, 
water or helium respectively). Therefore, the coolant will transmit thermal radiation over distances 
relevant to heat transfer in the core, and also absorb and emit it, depending on the salt composi- 
tion and state of contamination. The underlying theory and optical properties of relevant salts are 
described in Section 2.5. 


Flow, Convection and Conjugate Heat Transfer 


The behaviour of molten salts in terms of flow, pressure drop and convection is not particularly 
different to that of other more common Newtonian fluids. Salts do not exhibit complex rheological 
effects and do not have properties that are outside the range of others that are routinely analysed 
in engineering systems. Salts have relatively high Prandtl numbers (Pr), but other liquids (for ex- 
ample, oils and glycols) have substantially higher values, and water at room temperature is also a 
moderate Pr fluid. 


The behaviour in high Pr boundary layers is well understood, and well described by Kader (1981) 
and Kays (1994). The high Pr means that the thermal boundary layer will be smaller (or gradients 
higher) than the momentum boundary layer in laminar flows and steeper in the near-wall (sublayer) 
region of turbulent boundary layers (Figure 2.1). This is the opposite to the effect in fluids with low 
Pr, where the thermal boundary layer is much larger — this is discussed in detail in Volume 5 (Liquid 
Metal Thermal Hydraulics). The small thermal conductivity means that within the sublayer of the 
momentum boundary layer, turbulent energy transport plays a proportionally bigger role because 
of the smaller molecular transport. The effect of the high Pr needs to be taken into account in heat 
transfer correlations and in the near-wall treatment in Computational Fluid Dynamics (CFD), but 
the implications for modelling are not as systemic and important as they are for low Pr thermal 
hydraulics (especially for turbulent flows, because the effect is mainly confined to the near-wall 
region). 


The moderately high viscosity of salts and narrow flow paths in some core geometries means that 
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Figure 2.1: Temperature profiles for fully developed turbulent flow in a pipe (normalised 
by the difference between the wall, 7, and centreline, T-) at Re = 10* for a range of 
values of Pr, using the expression given by Kader (1981). The corresponding normalised 
velocity profile (U/U-) is shown for comparison, indicating close agreement in shape for 
Pr = 0.7, but steeper near-wall temperature gradients (smaller thermal boundary lay- 
ers) with higher Pr. Left: normalised radial position, right: comparing the same profiles, 
plotted logarithmically against y*. 


there will be more need to model low Re flows than for other reactor designs. For example, in chan- 
nel flow, there could be laminar-to-turbulent transition then re-laminarisation of flow downstream of 
an entrance or spacer. Alternatively, in a rod (or tube) bundle, the flow in the more open part of the 
cross-section could be weakly turbulent but the flow between the narrow rod-to-rod gaps could be 
laminar, resulting in complex and unstable interactions. 


Laminar-to-turbulent transition leads to a significant local increase in pressure drop and heat trans- 
fer, and this may manifest itself during a relatively minor transient or change in operating condition. 
More extensive transients could include reversal of the flow direction, for example in DRACS flow 
loops during a loss of forced flow event, and a transition to natural circulation. In these cases the 
effects of buoyancy on turbulence and velocity profiles in a channel, and the resulting impact on 
pressure drop and heat transfer, will need to be carefully accounted for. The guidance on these 
topics in Volume 3 (Natural Convection and Passive Cooling) is applicable to the flow of salts. In 
addition, there are flow paths in buoyancy driven reactor vessel cooling DHR systems that are 
external, so do not use salt as a working fluid, and for these the guidance in Volume 3 is also 
applicable. High Pr values, coupled with low Re flows and a substantial temperature dependence 
of viscosity, will mean that in some conditions there is a strong temperature and hence viscosity 
variation near to the wall of a channel where the temperature gradient is steep. This will need to be 
accounted for when applying heat transfer correlations. 


There are core designs, such as the MSFR (Figure 1.3), that have large open flow areas where 
power is generated, rather than more confined channel layouts. In these cases, there is an in- 
creased chance of complex effects caused by buoyant plumes or transiently occurring large turbu- 
lent structures, resulting in localised or rapidly changing variations in the temperature of the salt 
that impinges on a surface. This can have implications for structural integrity. 
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Some designs use pebble bed cores. While these have been studied in the context of HTGRs, 
the fluid properties of salts are likely to mean that the flow in the bed is in a different Re range for 
pressure drop and heat transfer predictions. An additional feature is that in FHRs, the pebbles are 
buoyant in the salt coolant (Andreades et a/., 2016); this is a feature that may require prediction in 
itself, and will also tend to dynamically redistribute the arrangement of pebbles in some regions of 
the core. This variability in possible geometric arrangements? will lead to localised uncertainty in 
flow and heat transfer performance in regions of the reactor. In addition, the localised variability in 
flow paths between pebbles will create hotter and cooler locations on individual pebble surfaces, 
and the fuel will need to be resilient to this. 


Liquid fuelled reactors feature a moving heat generating fluid and mobile heat generating species 
(fission products and DNPs). In the context of convective heat transfer, this is not an extensively 
studied flow configuration compared to the detailed assessment available in textbooks that has 
been carried out for convection in channels without internal heat generation. This is relevant to the 
applicability of convection correlations, because they implicitly contain information about the shape 
of the temperature distribution over the fluid cross-section; internal heat generation changes this 
shape. 


Conjugate Heat Transfer 


In common with other reactor types, MSRs require and will benefit from modelling Conjugate Heat 
Transfer (CHT). Including solid structures in an analysis provides the ability to move away from un- 
realistic, idealised boundary conditions and allows the stored energy in fuel and reactor structures 
to be incorporated. The guidance available in Volume 2 (Convection, Radiation and Conjugate 
Heat Transfer) is relevant and applicable to MSRs. 


A related phenomenon is the solidification of salt that could occur under either adverse conditions 
(such as an overcooling fault) or deliberately in long-term storage and material handling in filling 
or drain tanks. Incorporating realistic wall solids, potentially including nuclear heating, should be 
considered as the default option in these cases. For example, it was noted by Cartland Glover et al. 
(2019) that when CHT reactor walls were included in a MSFR reactor vessel simulation, compared 
to idealised boundary conditions, less freezing of the salt at a solid surface was predicted. 


Freezing is possible in primary circuit pipework and in heat exchangers, potentially preventing salt 
from flowing at all in some regions of a heat exchanger, or reducing the available cross-section 
for flow in channels? and creating a low conductivity solid wall layer. Another possibility for the 
modification of the heat transfer at the interface between a molten salt and a structural solid is 
the formation of high thermal resistance surface layers, especially where redox control is poor 
(Section 2.4). 


Considering solidification may also be necessary for pebbles or other solid fuel elements where the 
salt ‘wets’ the surface (depending on surface properties and the salt surface tension). Any adhering 


Aleatory uncertainty, see Volume 4. 

A range of more complex behaviours could also occur, such as freezing being inhibited in the centre of channels by the 
accelerated flow caused by narrowing of the cross-section (this may present issues for freeze plugs, where they cannot be 
fully frozen while the salt is flowing), or the formation of an oscillating flow where the salt repeatedly freezes due to localised 
cooling, but melts as a result of hotter temperatures that the flow reduction causes. 
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salt may freeze if removed to a lower temperature region (depending on the decay heat in the fuel), 
which will need to be considered in the design of refuelling processes and equipment. 


Thermophysical Properties 
There are two interpretations of the ‘variation’ of the thermophysical properties of a salt mixture: 


1. Variation of the properties of a salt of a given composition with temperature. 


2. The variation of the properties with the salt composition. 
This section describes the available data with respect to 1, and 2 is discussed in Section 2.4. 


There are a number of review papers and reports that collate and assess the available thermophys- 
ical properties for relevant salts. The incomplete and nuanced state of the available data means 
that obtaining and evaluating the applicable data in the following (and their references where nec- 
essary) is suggested: 


Cantor (1968) and Cantor (1973): Much of the data related to fluoride coolant and fuel salts orig- 
inates from these MSRE era reports. They contain important details on the experimental 
or estimation methods and uncertainties, qualifications and caveats that are not necessarily 
repeated or noted in more recent references that cite them. 


Janz et al. (1979) and Janz (1988): The first reference presents correlations for all thermophys- 
ical properties for 49 salt systems. The latter provides density, surface tension, electrical 
conductance and viscosity correlations (vs temperature) for hundreds of one and two com- 
ponent molten salt systems. While these are extensive resources, only a very small fraction 
of the salts are directly relevant to reactors, and this information (and that of many other ref- 
erences by Janz) is generally extracted and discussed in review articles below. However, the 
broader data relates to salt systems with constituents relevant to MSR coolants, and is used 
to help with property estimation methods for other salt compositions. 


Williams (2006a) and Williams (2006b): A wide ranging review of the characteristics and proper- 
ties of alkali fluorides, ZrF, salts, fluoroborate salts and chloride coolant salts, including the 
reasons to select them for a particular application. 


Romatoski and Hu (2017): A recent review of widely cited references which produces a set of 
recommended properties for LIF—BeF, (FLiBe), NaF—ZrF, (NaFZrF), LiF -NaF—KF (FLi- 
Nak), with representative uncertainties. No thermal radiation (optical) properties for these 
salts are presented and it is acknowledged that the dataset is incomplete, and that additional 
experiments to improve the confidence in the data and better quantify the uncertainties are 
necessary. For fluoride salts, this reference compares most of the others listed here. 


Khokhlov et al. (2009) and Benes and Konings (2012): A wide range of thermophysical, chemi- 
cal and thermodynamic properties for a range of fluoride coolant and fuel salts. Data from var- 
ious experimental and estimation sources are combined and discussed, with recommended 
values or correlations provided. 


Sohal et al. (2013) and Serrano-Lopez et al. (2013): Two reviews of a similar scope of the avail- 
able data for density, viscosity, specific heat capacity, and thermal conductivity properties for 
FLiBe, FLiNaK, LiF -NaF—BeF,, NaF—NaBF,, KF—ZrF4, KCI-MgClo, NaNO3;—KNO3 (So- 
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lar Salt) and NaNO,—NaNO,—KNOg (HITEC); the latter two are common in CSP systems. 
The need for more complete, high accuracy correlations with wider applicable temperature 
ranges (especially at higher temperatures) is noted. 


While the references above are able to provide property values or correlations for common salts, 
there would be considerable benefit in a consolidated database of improved properties (with more 
accuracy and over a wider range of temperatures) to support modelling and simulation. Jerden 
(2019) provides an interim report from a programme that is in the process of building a thermo- 
physical properties database. This report reviews availability of a range of salts and concludes: 


...some key properties for MSR-relevant salts exist in the open literature. However, 
the existing literature data sets are sparse and incomplete for the actinide-bearing 
salts relevant to modern MSR applications. It is concluded that there remains a need 
for an experimental program focused on collecting data needed to address the design 
and licensing needs of current MSR concept developers. 

... There is also a need to continue to strengthen the technical link between the MSR 
modeling and simulation work and the experimental thermophysical properties mea- 
surement work. For example, MSR models could provide sensitivity runs to rank the 
relative importance of particular salt systems and thermophysical properties so that 
the experimental work can be prioritized and more directly linked to the modeling and 
simulation work needs. 


Example Salt Properties 


Salts that are likely to be used in reactor applications can be characterised, in comparison to water, 
as having similarly high volumetric heat capacities (pcp, one of the few coolants where this is the 
case), similar thermal conductivity, and substantially higher viscosity, giving higher Prandtl numbers 
(Pr = 10 to 30). These properties are compared to water (and helium where suitable) for a number 
of coolant and fuel salts in Figure 2.2 over the approximate temperature range in which they would 
be employed‘. The melting points, Tjy, of the salts are also shown in the legend. 


In addition to the basic thermophysical properties (p, w, Cp, k, pc, and Pr), the temperature depen- 
dence (at constant pressure) of density (volumetric thermal expansion, 6) is plotted: 


This is 6 that appears in the Grashof and Rayleigh numbers, and so is relevant to buoyancy driven 
flow and buoyancy effects in heat transfer. The values for B in the salts shown are similar to each 
other, and are substantially lower than for water or helium. 


The normalised viscosity gradient 
1 Ou 


p OT 
is also shown, which for salts is similar to that of water (at reactor temperature). This is relevant 
because, in channels where there is a substantial temperature variation over the cross-section (due 
to internal heat generation and/or convection transfer to the channel walls and as a result of the 


4 For water this is from 25°C to close to saturation temperature at 155 bar (typical PWR primary circuit pressure). 
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Figure 2.2: Thermophysical properties of example salts, water and helium. 

Data sources: FLiBe, NaFZrF, FLiNaK: Romatoski and Hu (2017) and Bene and Konings (2012); 
LiF—ThF: Rouch et al. (2014); LIF -BeF,—ThF,—UF,: Cantor (1968); 

NaNO3—KNOs;: Serrano-Lopez et al. (2013) and Jin et al. (2016); 

Water: IAPWS (Wagner et al., 2000); He: NIST Chemistry WebBook (webbook.nist.gov). 


17 of 89 


Volume 6 


Technical Context o; 


Volume 6 


high Pr), then there will be a corresponding variation in viscosity, which may need to be accounted 
for in a convection heat transfer correlation (Section 3.1). 


A similar set of figures and more detailed comparison to a wider range of coolants, including su- 
percritical water, COz, and liquid metals (but only one salt, FLiNaK) can be found in Appendix A.2 
of Pioro (2016). 


There are similarities in the responses to temperature of the example salt properties which are 
described below. The behaviour and functional forms indicate the variation of a property that should 
be expected with temperature. These descriptions are intended to assist in interpreting Figure 2.2 
as well as the property estimation methods described in Section 2.3.2. They are also the forms 
of the equations that are typically found in the evolving databases and literature regarding salt 
properties described above. 


Melting point: Fluoride salts have melting points in the region of 500°C; the nitrate solar salt is 
significantly lower at 222°C. 


Density: The density of each salt has a linear decrease with temperature: 
p=A-BT 


Viscosity: Molten salts exhibit Newtonian viscosity behaviour (Holcomb and Cetiner, 2010), and 
the typical exponential decrease with reciprocal temperature (which can be referred to as the 
Andrade equation): 

i AeB/T 

Specific heat capacity: Thermochemistry data in the NIST Chemistry WebBook’ is provided as 
functions for enthalpy and entropy at standard pressure as a function of temperature. These 
can also be formulated to give c,: 


G =A+BT+CT?+ DT? +E/T* 


Coefficients for this equation are available for pure substances, including salt constituents of 
interest, but not for mixtures. Because the measured accuracy of this property will be limited, 
it is unlikely that more than the linear terms (A and B) will be used, although the database 
under construction at Argonne National Laboratory (ANL) (Jerden, 2019) does also have 
provision for the terms in T? (A, B, C and E). 


Thermal Conductivity: The hardest property to measure for molten salts (Cooke, 1973), which 
led to some confusion and errors in heat transfer calculations in the MSRE programme. 
Where the temperature dependence is presented, it is linear with temperature: 


k =A+ BT 


5 webbook.nist.gov, originating from Chase (1998) 
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Property Estimation Methods 


The uncertainties in the salt properties shown in Figure 2.2 are large compared to those for wa- 
ter, despite these being salt mixtures (in terms of their species and concentration ratios) that are 
amongst the most well-studied. In some cases, the values of c, and k are plotted as being con- 
stant with temperature; only single values are presented in the references, largely because the 
uncertainty in the value is of a comparable magnitude to the expected variation with temperature. 


In reactor component development and testing, it is likely that alternative/simulant salts will be 
used, especially where fissile materials are to be dissolved in the final salt. The composition of 
a salt may also vary with irradiation or operational adjustments. Therefore, it is not expected that 
experimental data for the properties of all compositions of a salt that may be encountered will be 
available. 


For salts without reliably measured data, the issue becomes how accurately it can be determined, 
either by derivation from first principles, or by estimation from correlations derived from measured 
data for similar salts (or the same salt constituents at a different concentration ratio). Directly mea- 
suring the properties is relatively difficult because the measurements need to be made at high 
temperatures and there are challenges in purifying and handling materials correctly to obtain ac- 
curate measurements. Without measures to control contamination and corrosion, the data is easily 
corruptible. This requirement for purity does, however, indicate that in reactor conditions, where 
ideal laboratory conditions are not likely to be achievable (nor necessary), then there will be a level 
of uncertainty or variation compared to results for pure substances. 


The methods available for estimating values for the five thermophysical properties discussed above 
(Section 2.3.1) are described below. Guidance is available on the evaluation of density, viscosity, 
heat capacity, and thermal conductivity for various molten salt mixtures Na, Be, K and Li fluorides 
with Zr, Th, U, and Pu fluoride components. Less information is available for chloride salts, al- 
though in many cases the methods described should be equally applicable, and can be compared 
to the data that does exist (such as Bystrai et a/., 1974 and Katyshev ef al., 1982). The details of 
the methods are not reproduced here because of the level of complexity and nature of the quali- 
fications, caveats and notes on accuracy and applicability. The noted references should be used 
directly. 


Melting point: According to Holcomb and Cetiner (2010), an extensive database of phase dia- 
grams exists for salt systems of all types, and there is no need to pursue estimation tech- 
niques. The majority of the binary phase diagrams of interest have been measured, as well 
as many of the ternary systems. Sources for this data are described in Section 2.4.3. 


Density: According to Williams (2006a), density is among the most straightforward of properties 
to measure and is one of the most readily estimated for new compositions. The estimation 
method is a mole fraction weighted sum of the molar volume of the constituent salts in a mix- 
ture, which themselves are known from experiments. Given that the variation of density with 
temperature is linear, if molar volumes of the constituents are known at two temperatures, 
then the temperature variation of the mixture can also be calculated. This method is also 
applied by Khokhlov et al. (2009) and Benes and Konings (2012), where examples are given 
and the accuracy discussed. 
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Viscosity: Viscosity is the property that varies most with temperature and salt composition. There 
are no simple methods for estimating it from first principles, although many systems have 
been studied experimentally. Determining the viscosity (and other properties) of an arbitrary 
salt mixture by molecular simulations is, however, a developing method, and will be discussed 
in Section 4.1. Khokhlov et a/. (2009) describes a method for combining experimental data 
for the molar viscosity of salt components and (preferably) binary mixtures based on their 
molar proportions to estimate the properties of a ternary mixture. For example in Khokhlov 
et al. (2009), the viscosity of two compositions of LiF-BeF,—ThF, are estimated from the 
viscosity of the LiF—BeF, and LiF—ThF, binary mixtures. Again, Benes and Konings (2012) 
give examples and comments on accuracy. 


Specific heat capacity: The molar mixing methods that are suitable for density and viscosity can 
also be applied for c, where the data of the constituents is known. Because there is a lack of 
suitable experimental data, it cannot be applied for the heat capacity of mixtures containing 
Zr, Th, U, and Pu fluoride. However, an empirical linear correlation between experimental 
values for cp (J/kg K) for 30 molten fluoride salts (including those containing actindes) and 
inverse molar mass (M, kg/kmol) was found by Khokhlov et al. (2009): 


Cp = 291.6 + 8.02 x 10*/M 


Bene§ and Konings (2012) note the following about this method: 
Data on experimental heat capacity of molten fluoride systems containing ac- 
tinide fluorides are generally lacking. [Khokhlov et al. (2009)] derived an estimation 
equation that is based only on the reverse molar mass of the salt mixture. However, 
this method must be considered only as an approximation as it sometimes gives 
relatively large deviations to the experimentally determined heat capacity (e.g., 
pure LiF or NaF). Thus, more measurements are required in order to describe the 
heat capacity behavior of the liquid fluoride solutions more precisely. 
An alternative is the Dulong and Petit method, which assumes a constant contribution to 
Cp per mole of each atom in the mixture. It is described and compared to experiment and 
the above expression of Khokhlov et al. (2009) by Williams (2006a) and Romatoski and Hu 
(2017). Computational thermodynamics methods (Section 2.4.5) are also able to provide c, 
assessments. 

Thermal Conductivity: The lack of reliable sources of experimental data for k in salts with rele- 
vant compositions makes the estimation of more complex mixtures by molar mixing not pos- 
sible. However, Khokhlov et al. (2009) does provide an empirical expression for k (W/m k) 
based only on temperature (K) and molecular mass (kg/kmol): 


k =—0.34+0.5 x 10-37 + 32/M 


Bene§ and Konings (2012) also provide a discussion of this method: 
Experimental data for the thermal conductivity of molten fluoride systems are gen- 
erally lacking. Although [Khokhlov et al. (2009)] recently derived an estimation 
method for this quantity on the basis of the available data for the molten chlo- 
rides, bromides, and iodides, more measurements are needed in order to justify 
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his approach. 
An alternative is the Rao-Turnbull correlation that depends on molecular weight, salt density 
and melting temperature, and again is discussed and compared by Williams (2006a) and 
Romatoski and Hu (2017). 


Salt Chemistry, Thermodynamics, Dissolved Gases and 
Nuclear Effects 


The scope of guidance on methodology described in Section 3 is limited to predictions of flow and 
heat transfer under single-phase (liquid) conditions, and does not provide details on the methods 
described in this section. However, many of the phenomena of interest in analysing MSRs are not 
thermal hydraulic per se (if that is defined as pertaining to the temperature, pressure and velocity 
fields of a fluid), but they either: 


* depend on temperature, and so require accurate heat transfer assessment, or 


* are a function of salt composition, which will define (or introduce uncertainty into) the ther- 
mophysical properties that a thermal hydraulic calculation or simulation will use, or 


* involve species transport and potentially reactions within the primary, secondary or DHR 
circuits. 


Therefore, it is necessary to be aware of how any thermal hydraulic analysis fits within or interfaces 
to this broader context. In many cases, these broader issues could be implemented as a sub-model 
within the structure of a thermal hydraulic tool, or in another tool that is closely coupled to it. 


Salt Chemistry 


For liquid fuelled reactors, the composition and chemistry of the salt is the issue around which most 
other considerations revolve: 


Fuel chemistry is an area of vital importance to a fluid-fuel reactor comparable to fuel 
structure, cladding integrity, and coolant stability in a solid-fuel reactor. 


Haubenreich and Engel, 1970 


Chemical considerations are key to the selection of a salt to use as a fuel or coolant, and monitoring 
and controlling its chemical state will comprise a large portion of the instrumentation, operation and 
maintenance challenge of any MSR. One of the principal considerations is material compatibility; 
graphite is generally not reactive with the salts chosen when it is used as a fuel component or 
moderator structure, so the reactivity of the salt with structural metals (corrosion) is of greatest 
concern. 


The alloy Hastelloy N was developed by the MSRE programme. Using this material, it was possible 
to control corrosion sufficiently well to have confidence that it would not be a significant technologi- 
cal barrier to reactor development (Grimes, 1970, Haubenreich and Engel, 1970). These materials 
and the understanding of their corrosion behaviour have been the subject of substantial ongoing 
research (Delpech et al., 2010, Zheng and Sridharan, 2018). It is accepted that it is necessary 
to exclude oxygen and moisture from a salt and that active ‘redox control’ should be applied to 
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minimise corrosion (Zhang et a/., 2018). The redox potential of a salt measures its tendency to ac- 
quire electrons from or lose electrons to structural metals and be reduced or oxidised respectively. 
Holcomb and Cetiner (2010) discuss the reasons why corrosion control is very important, and an 
overview of corrosion in the context of broader chemical considerations can be found in Ottewitte 
(1992) and Fredrickson et al. (2018). 


A corollary to corrosion and dissolution of structural materials is the deposition process, where 
noble-metal fission products can ‘plate out’ on surfaces (Haubenreich and Engel, 1970). This ef- 
fect and the solubility and oxidation state of other salt constituents depends on the redox conditions 
(which are affected by the presence of fission products). Fission products can detrimentally alter 
the neutronic performance of a core, and so the design and operation of reactor chemistry may 
want to enable the extraction of certain elements from the salt during operation. Rigorous predic- 
tion of these processes requires transport modelling of their production rates (or thermodynamic 
equilibrium concentrations) and of surface deposition reaction rates, all of which can have a depen- 
dence on temperature. This could be achieved in coupled, multiphysics simulations underpinned 
by thermal hydraulics, although bounding and exclusion of the effects by simpler methods may be 
preferable (at least during reactor concept design). 


For statically fuelled reactors, a subset of the salt chemistry issues are equally important, but other 
aspects are less severe because of the comparative stability of the composition of the fluid and its 
significantly reduced radioactivity and heat generation. 


From a thermal hydraulic perspective, the most directly relevant phenomena are: 


* The development of high thermal resistance layers on salt-metal interfaces where redox con- 
trol has not been implemented and deposition occurs, which will alter convective heat transfer 
performance. 


¢ Dissolution of constituents of structural alloys (such as chromium) or the presence of fission 
products can significantly alter the optical properties of a salt, changing the radiative heat 
transfer behaviour. 


Nuclear Effects 


Nuclear effects mainly impact on thermal hydraulics modelling in the form of heat sources in solid or 
liquid fuel regions, and nuclear heating from neutron moderation and absorption as well as gamma 
absorption in the salt, moderator and structure. A PIRT assessment of the neutronic modelling 
phenomena and needs can be found in Rahnema et al. (2019) (solid fuel) and in Diamond et al. 
(2018) (liquid fuel). 


In all types of MSR there is an important nuclear reactivity feedback from increases in fuel temper- 
ature as a result of Doppler broadening. For liquid fuelled reactors, as described by Krepel et al. 
(2014), there is a further strong negative feedback from fuel density reduction, and in a thermal 
spectrum reactor, also potentially from moderator expansion (reducing the liquid space available in 
the core). This strong interdependency requires the physics of neutronics and thermal hydraulics 
to be coupled for operational and safety analysis in liquid fuelled reactors. This needs to account 
for the motion of DNPs, and the heat generation within the fluid from the decay of fission products, 
which will produce heat outside the core region (for example, in heat exchangers). 
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The presence of fission products and transmuted elements in the salt modifies its composition, 
which itself modifies the neutronics behaviour. The new species, and changes in valence state 
caused by their production, impacts on the salt chemistry, and the solubility and retention of 
gaseous fission products or the deposition of other solids, including fissile constituents. The im- 
plications for transport, retention and behaviour of these radiotoxic elements therefore is linked to 
Sections 2.4.1 above and 2.4.3 below. 


Melting and Solidification: Salt Mixture Phase Equilibria 


It is unlikely that a single composition salt would be used in a reactor application; a mixture of salts 
is used to reduce the melting point. An example of the variation in melting point with composition 
is shown in Figure 2.3, the phase diagram of the binary mixture LiF—BeF. (a common choice of 
salt, Table 1.2). Starting at the left, with pure LiF, the decrease in the liquidus® temperature (melting 
point) with increasing BeF, concentration can be seen. It reaches a peritectic equilibrium point at 
a 33.5 mole % BeFs, 454°C, then a eutectic equilibrium point at a 52 mole % BeFo, 355°C, before 
increasing with increasing BeF». Note that more recent evaluations of this salt have produced 
slightly different equilibrium points for this salt system compared to the 1950s ORNL data shown 
in Figure 2.3, hence the melting point of 456°C shown in Figure 2.2. 


To have the minimum melting point, the eutectic composition would be chosen, but this salt is 
a good example of why this is not always the case. For example, one of the constituents may be 
significantly more expansive than the others, or there are implications for properties such as vapour 
pressure or viscosity. Williams and Clarno (2008) notes (citing data in Williams, 2006a) that having 
BeF, as a constituent in salt mixture can cause a large increase in viscosity: 


In general, the other salt properties do not place severe limitations on the choice of 
coolant composition, and they are usually considered in combination with respect to 
their influence on heat transfer. The one exception to this concerns the viscosity of 
BeF.-containing salts. Salts rich in BeF, form extremely viscous mixtures that com- 
pletely restrict their use to compositions of BeF, < 45%. Most other salts do not pos- 
sess this property of forming highly viscous mixtures. 


Similarly, because of substantial increases in vapour pressure: 


... the range of useful compositions is usually restricted to 40 to 45% of BeF» or ZrF 4, 
depending upon coolant 


hence the 66.5—33.5 mole % LiF—BeF, peritectic composition is favoured, despite the higher melt- 
ing point compared to the eutectic composition. One implication of this is that, in low temperature 
conditions, at the onset of solidification (dropping to the liquidus), one salt component will preferen- 
tially crystallise and the remaining liquid salt composition will be altered. For simple closed systems 
under equilibrium conditions (the type of system used to generate phase diagrams) this process 
will be reversible and the total composition is not altered. However, in a complex flow process in, 
for example, a heat exchanger, this may lead to preferential deposition on surfaces of one species 
and localised composition, and hence property and potentially reactivity changes. 

The liquidus is the line separating the area of the phase diagram that is all liquid from that where there is a mixture of 
liquid of one component as well crystals of another. The solidus is the line separating the area of the phase diagram that 


is all solid from that which is a mixture of liquid and crystals. For a more detailed introduction to the terminology, see 
www.tulane.edu/~sanelson/eens21 1/2compphasdiag.html. 
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Figure 2.3: Phase diagram for LiIFp—BeF, (FLiBe) (Weaver et a/., 1961). 
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Figure 2.4: Phase diagram for LiF—NaF—KF (FLiNak) (Thoma, 1959). 
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Phase diagrams can become very complex with more than two components and multiple equilib- 
rium points; a simple ternary system is shown in Figure 2.4, which is a projection of the liquidus 
surface against the concentration of the three components. It has a single eutectic point at 46.5— 
11.5-42 mole % LiIF—NaF—KkF, 454°C. 


Phase diagrams can be found for a range of salts considered for reactor application in Thoma 
(1959), Weaver et a/. (1961), Sangster and Pelton (1987) and Benes and Konings (2012). Ricci 
(1958) provides guidance on their interpretation. Data on the eutectic composition and melting 
point of more than 6000 salt mixtures is given by Janz et al. (1978). 


Dissolved Gases, Boiling and Vapour Pressure 


One of the reasons to choose a molten salt as a reactor coolant is their high boiling point’. In the 
candidate salts for reactors, the boiling points are very high (>1400°C, Williams, 2006a, Benes 
and Konings, 2012) and there is therefore expected to be large thermal margins to boiling, even for 
onerous postulated accidents. 


While boiling is not a significant concern, the vapour pressure of salts is of interest in designing 
reactors that are likely to have pool-like layouts (with a free surface) in either their core or other 
system components. This will require an inert cover gas, and gaseous fission products and tritium 
are expected to require removal, as described by Williams (2006a): 


Even in a low-pressure system, the magnitude and nature of vapor produced from the 
salt needs to be evaluated. Experience with the ARE and the MSRE shows that very 
low salt vapor pressures (<1mm Hg) simplify the off-gas system design and that cer- 
tain vapor species can present problems. 

In any high-temperature salt system, a purged cover gas will be necessary. The trans- 
port of significant amounts of salt vapor in this cover gas system can cause problems. 
In the operation of the ARE, it was found that the vapor over the ARE salt (63%NaF— 
41%ZrF4-6%UF 4) was nearly pure ZrF 4. Because ZrF, sublimes rather than boils, ZrF4 
‘snow’ was found in the exhaust piping. 


The composition of the salt mixture affects the vapour pressure considerably. Not all components 
of a mixture will have the same vapour pressure, and the fact that they are in mixtures alters 
their vapour pressures (their ‘native’ volatility is suppressed). It may also be the case that some 
dissolved species that arise as fission products will be volatile at normal operating conditions, while 
the underlying carrier salt is not. 


Correlations for the variation of vapour pressure with temperature for some common MSR mixtures 
are given by Bene§ and Konings (2012); for example, for LIF -BeF, (66-34 mole %) in the range of 
823 to 1473 K the vapour pressure, Pyap (in N/m?) as a function of temperature, T (in K) is 


P... = 10 (11-914—-13003/T) 
vap — 


giving, Pyap = 6.7 N/m? at T = 1173 K (900°C). 


7 The temperature at which the partial pressure of vapour salt species above a liquid surface is 1 atm. 
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The retention and/or controlled release of gaseous fission products (such as xenon and krypton) 
will need to be addressed in any design, and will depend on temperature. Some data exists for 
solubility (e.g. IUPAC, 1991). It was shown in the MSRE that xenon could be effectively ‘stripped’ 
(removed) from a fuel salt into an off-gas system during operation by sparging (bubbling) an inert 
gas into it, although this did have the effect of introducing a ‘mist’ of fine salt droplets into the gas 
space and allowing gas bubbles to circulate in the core, introducing perturbations in reactivity due 
to voids (Haubenreich and Engel, 1970, Engel et a/., 1970). In a system where entrained gases 
are encountered, knowledge of the surface tension of the salt will be required, and is available for 
a number of salts (Desyatnik et a/., 1975, Janz, 1988). 


Another dissolved gas that requires consideration is tritium, which is produced from several neu- 
tronic reactions with Li, Be, F, in the fuel and coolant salt, and by ternary fission. At high tempera- 
tures, hydrogen and (hence tritium) can diffuse through structural metals. It can therefore escape 
from the primary circuit, principally through heat exchangers because they have a deliberately large 
surface area and thin walls as part of their design. This means that tritium will migrate to secondary 
or DHR circuits. There are methods available to control the concentration and migration of tritium, 
and they can interact with the redox control of the salt. Gas sparging is able to remove tritium from 
molten salts in the same way as with Xe and Kr, and tritium can be captured onto materials such 
as carbon. A review on the current state of knowledge and development is given by Forsberg et al. 
(2017). 


Computational Thermodynamics 


Experimentally determining the properties of mixtures of salts is not always practical or economic. 
To address this, computational thermodynamic modelling of phase diagrams for multi-component 
mixtures (a semi-empirical calculation method referred to as CALculation of PHAse Diagrams 
(CALPHAD), Saunders and Miodownik, 1998, Pelton, 2019) is used. This takes information that 
is known about the pure constituents of the mixture® to predict the phase diagram of complex 
mixtures. This allows the assessment of melting behaviour as well as the possibility for the precipi- 
tation of solid phases (including noble-metals or fission products). The solubility of dissolved gases, 
vapour pressure of the salt, redox potentials, density and specific heat capacity are also able to be 
predicted as a function of salt composition and temperature (Besmann and Schorne-Pinto, 2021). 


This prediction is made by using a numerical Gibbs Energy Minimisation (GEM) algorithm to deter- 
mine the interfaces between phases; Gibbs energy is related to enthalpy, temperature and entropy 
by 

G=H-TS (org =h-—Tsas‘specific’ quantities) 


For a condensed (liquid or solid) phase of a pure constituent (Smith et a/., 2020) 


T Te 


é(T)dT — Ts* rf dT 
0 


B(7) = ark? + [ 
To 


where g°(T) is the molar Gibbs free energy of the material in the form of a potential, with units of 
J/mol. The enthalpy of formation, Arh®, is the change in the sum of the internal energy and flow 


Pure constituents are known as ‘end-members’ — they are the ends/corners of phase diagrams, so LiF and BeF, are the 
end members on Figure 2.3. 
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work (PV) resulting from the formation of 1 mol of the substance from its constituent elements at 
the standard state, which is a pressure of 1 bar and a temperature of To = 298.15 K (25°C, for more 
detail see Rogers and Mayhew, 1992). The absolute entropy at the standard state is represented 
by §* (units JK~! mol~}), which is the change in the entropy that results from increasing the 
temperature of a substance from 0K to 298.15 K. Here ¢, is the molar specific heat capacity at 
constant pressure (also units JK~! mol~+). This description only accounts for dependence on 
temperature (because it is typically dominant), but it is possible to include the effect of pressure 
dependence on the Gibbs energy. Data for A-h®, §° and € are available from the NIST-JANAF 
Thermochemical Tables (Chase, 1998)° for a very large range of substances. 


The Gibbs energy function that results from a specific mixture of such components and the ther- 
modyanmics of their interactions determines its response as a function of temperature, and the 
equilibrium state is found by minimising it (maximising entropy). Commercial GEM tools such as 
FactSage, Thermo-Calc or open source tools such as OpenCalphad or Thermochimica are avail- 
able to perform this modelling, and the availability of databases for the necessary inputs are a 
key consideration. The application of CALPHAD methods is common in the assessment of metal 
alloys (Kattner, 1997), where there have been extensive studies, and the necessary models and 
data exist. The nature of salts as ionic liquids, that display local structure (short range order) mean 
that the models and data are less extensive and mature. 


Computational thermodynamics is one part of a set of tools and specialist scientific disciplines 
needed to understand the detail and reasons underlying the behaviour and properties of salts. It 
is necessary to combine available data and comparisons from all sources and aggregate them 
to construct the best model possible. To provide support to modelling and simulation, ORNL are 
creating a Molten Salt Thermodynamic Database (McMurray et al., 2019), and it is acknowledged 
that more measured data is needed to be able to reliably make assessments (McMurray et al., 
2021). A similar project is underway at the European Commission Joint Research Centre (JRC) 
(for example, Schacherl et a/., 2021 and Bene§ et al., 2021). While some data sources are experi- 
mental, there is also an increasing use of molecular modelling to supply inputs, and to complement 
CALPHAD type calculations (Smith et al., 2020); this is discussed in more detail in Section 4.1. 


Radiative Heat Transfer 


Radiative heat transfer is a complex and nuanced discipline of physics in its own right, and it is 
generally only applied in a simplified and approximate manner in mechanical engineering anal- 
ysis. In systems containing water or liquid metal, it is not considered because thermal radiation 
does not penetrate any appreciable distance before it is absorbed (although radiative heat trans- 
fer is highly relevant in the inert gas space above, for example, a liquid metal pool). In systems 
containing non-participating media (one that neither absorbs, emits or scatters radiation at the In- 
frared (IR) wavelengths of interest, such as air at moderate temperatures), there are established 
methods for predicting radiative heat transfer between surfaces, where simplifications to ignore 
the variation of surface properties with temperature and wavelength are usually applied. Radiation 
plays an important role in heat transfer in combustion systems, where water vapour, CO, and soot 
are present at high temperatures and in significant concentrations, and emit, absorb and scatter 


9 and from webbook.nist.gov 
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radiation significantly. In a nuclear thermal hydraulics context, for heat transfer over moderate dis- 
tances, pressurised CO, (as present in AGRs), steam or gas mixtures with significant water vapour 
concentrations may need to be considered as a participating medium. 


From this starting point, the following questions are relevant to identify what is different about 
thermal radiation in molten salts, and how well it is understood: 


1. Are salts transparent? Do they participate in radiation? Do they scatter radiation, or just emit 
and absorb? 


2. What components, flows, length scales and temperatures require radiation to be considered? 
3. How does the emission and absorption of a salt vary with temperature and wavelength? 


4. How do specific salts compare to each other and to other materials in terms of absorption 
and refractive index? 


5. How complete and accurate is the available material property data for relevant salts? 


6. What effect does the surface have (emissivity, changes with ageing, corrosion, contamina- 
tion)? 

7. What is the effect of additives like fissile materials or fission products (such as actinides, lan- 
thanides or transition metals) or contaminants from impurities from the corrosion of structural 
materials? 


8. What additional considerations are necessary beyond (or approximations necessary to use) 
conventional heat transfer based on the familiar q = eoT* relation (as in Volume 2, Sec- 
tion 2.1.3; € is the surface emissivity’? and o is the Stefan Boltzmann constant). 


In general, simplified (and where possible bounding) models for accounting for radiation will be 
preferred. When these are implemented in a CFD code, the solution of the Radiative Transfer 
Equation (RTE) will already be implemented by the software, but the question will be, ‘what are 
the appropriate models, settings and properties to apply?’ For simpler calculations and system 
code models, similar choices will need to be applied, but in a form appropriate for the model. 
The approaches to implementing these calculations will be described in Section 3; the underlying 
theory and methods for developing the appropriate material properties are discussed below. 


Spectral Variation and Absorption 


To understand how to interpret the optical properties of a salt to calculate effective coefficients to 
be used in radiative heat transfer predictions, it is necessary to understand both: 


* how the spectrum of emitted and absorbed wavelengths varies with temperature, and 
« how this is combined with the optical properties of a salt that themselves vary with wave- 
length. 


The blackbody emissive power from a (black) surface surrounded by a transparent medium of 
10 Modest (2013) states that strictly the term ‘emittance’ should be used for real surfaces with ‘emissivity’ reserved for pure, 
perfectly smooth surfaces. The same applies to absorptance vs absorptivity, reflectance(-ivity), transmittance(-ivity). How- 


ever, there is not a consistent application of terminology, and the -ivity convention familiar from engineering references such 
as Incropera et al. (2011) will be applied. 
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constant refractive index, n, is a function of temperature and wavelength, as given by Planck’s law 


2mhca 


Ep = n2d5(eheo/mAkeT — 1) 


where k, is Boltzmann’s constant, h is Planck’s constant and cp is the speed of light in a vacuum 
(Modest, 2013). It has the units of W/m? um, so is a spectrally resolved energy flux emitted from 
a surface (subscript \ denotes a wavelength varying (spectral) quantity). The total emissive power 
of a blackbody is 


E, = i. E,,d\ = n’°oT* 
0 


where it can be seen how the Stefan-Boltzmann constant arises from a combination of the funda- 
mental physical constants found in Planck’s law: 


2n° kh 
o = —~3 
15f2c2 


A related quantity that will be important in the RTE, because it accounts for the directionality of a 
radiation field, is the radiative intensity, I, which is the energy emitted from a surface in a particular 
direction (so the units are W/m? sr um, where sr is steradian, solid angle). Integrating the blackbody 
radiative intensity leaving a blackbody, I, over the solid angles in a hemisphere gives 


Ipy = Epy/0 


and the total intensity is the integral of the spectral intensity over all wavelengths: 


ii Ib dX 
0 


For real surfaces, their emissivity, absorptivity and reflectivity depend on surface temperature, and 
the direction and wavelength of the radiation as well as the surface finish or condition. It is com- 
mon and preferable to average over all directions and wavelengths to obtain a single value for the 
hemispherical emissivity. Another common approximation, that leads to a single value for a surface 
property, is that surfaces are diffuse (there is no preferred direction) and ‘gray’, where there is no 
variation with wavelength over the region of spectrum where E,) is large, and that the property 
does not depend on temperature. A consequence of applying the assumption of local thermody- 
namic equilibrium, stated by Kirchoff’s law, is that the emissivity and absorptivity of a surface (either 
as single values, or resolved by temperature, directional and spectral dependence) are the same, 
so € = a can be applied. This also applies in the volume of any participating medium. From this, 
if scattering is ignored, the quasi-steady form of the RTE for the radiative intensity in an absorbing 
and emitting medium"', along a given direction for distance s, at a given wavelength is 
dIy 


we Ky(Ip, — Ty) 


where x is the absorption coefficient (with units 1/m). Integrating with respect to the direction ds 


"1 It is not possible for a medium to absorb radiation but not emit. 
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from 0 to s for a isothermal layer gives 


I,(s) = L(O)e™ + Ioa(1 — 7) 


Ss 
TT -/ K ds 
0 


and is known as the optical thickness (for a spatially uniform value of k, then 7, = kys). k is 


where 


related to k, the imaginary part of the complex index of refraction, m 
m=n-—ik 


by (Modest, 2013, Equation 2.43) 
kK = 4anok 


where 7 is the wavenumber in a vacuum (7 = 1/Ao). 


Considering only absorption, and assuming no variation with wavelength, the more familiar quan- 
tity of absorptivity, a, from Beer’s law is found (omitting for simplicity, but this equally applies 
spectrally): vila 
= 4S 
_ O18) = _ 


Similarly, considering only emission (I(0) = 0), emissivity is defined as the ratio of the energy 


emitted to emission by a blackbody at the same temperature 


e=——=1-e” 


which, as for surfaces, is the same as the expression for absorptivity. 


Solving the RTE in general is difficult, partly because it needs to be integrated over all angular 
directions, in three dimensions of space, and over all wavelengths (or alternatively be framed as 
gray — without wavelength dependence in properties). It requires the temperature of the medium at 
all points to evaluate the blackbody function, I,, which requires the thermal radiation to be coupled 
to and solved with the overall energy equation for the medium. Solving the RTE can be achieved in 
one dimensional geometries with the gray assumption, or included in combination with conduction 
and convection, by approximations requiring a domain to be an optically thick (rt >> 1) or thin 
(7 < 1) geometry. It can be solved numerically in realistic 3D geometries by a number of methods 
included in CFD codes (Section 3.3). 


2.5.2 Thermal Radiation Material Properties 


The IR absorption spectrum of liquid water is shown in Figure 2.5, and is compared to the black- 
body emission spectra (Planck’s law) at two high-temperature reactor relevant temperatures (500°C, 
1000°C), and the solar spectrum. The absorption coefficient of water increases quickly after the 
end of the visible spectrum, and is large over the reactor temperature wavelength ranges (the tem- 
peratures relevant to water cooled reactors are even lower, hence further into the highly absorbing 
range). 
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For molten salts, the question is ‘how does the absorption spectrum compare, how much does this 
depend on salt composition and temperature, and how much data is known?’ Significant research 
has been carried out into a range of single composition crystals that are constituents of reactor 
relevant salts because they are used as optical components for high-powered laser systems with 
wavelengths in the near-IR". 


The Lorentz model is the relevant theory for the interaction between electromagnetic waves and 
ionic crystals. It predicts absorption in the Ultraviolet (UV) region caused by electronic excitations 
by short wavelength photons (referred to as the ‘Urbach tail’) and strong absorption at long wave- 
lengths due to photon excitations of lattice vibrations (phonons) which begins in the mid-infrared 
(also called ‘Reststrahlen bands’). In between these regions, the material is transparent, with the 
remaining absorption attributable to many factors, notably crystal vacancies and dislocations (in 
solid crystals), impurities and surface contamination. 


Li (1980) produced a correlation and coefficients for the phonon absorption region of alkali halide 
salts, including LiF and NaF and NaCl'’. A schematic representation of the absorption spectrum 
of an alkali halide crystal based on the Lorentz model is shown in Figure 2.6. The wavelengths 
where the UV and IR bands (steep increases) occur depend on the material, and the absorption 
level in the impurity dominated ‘floor’ will depend on the type and concentration of impurities. The 
correlation of Li (1980) was used, along with other experimental data sources by Derdeyn et al. 
(2018) and Chaleff et al. (2016), who produced estimates of the absorption spectrum of FLiBe and 
FLiNaK respectively. These estimates are shown in Figure 2.7, along with a comparison to water 
and to LiF and NaCl for reference. It can be seen that the wavelengths where absorption starts 
to increase significantly includes a substantial portion of the wavelength spectra at the reactor 
relevant temperatures (particularly for FLiBe). 


To add complexity, as the salt or surface temperatures increase, the emitted wavelengths move 
towards the low absorption region, and so a salt will become more transparent to thermal radiation 
with increasing temperature. This may have an effect on the progress of a reactor event, where 
increased temperature will cause less of the emitted radiation to deposit within the coolant and 
will increase the exchange between solid surfaces. Because a large proportion of the heat trans- 
fer is in the impurity dominated ‘transparent’ wavelength region, the effect of impurities (such as 
chromium from structural materials) may significantly increase the absorption in this region, affect- 
ing the overall absorption (Chaleff et a/., 2016, Chaleff et a/., 2018). For salts that are used in CSP 
applications (including solar nitrate salts and (K—Na)Cl), their absorption spectrum is such that 
they absorb more strongly in the visible wavelength region (Tetreault-Friend et al/., 2017) and can 
have their absorption deliberately and substantially modified, for example by the inclusion of small 
concentrations of other metal chlorides (Myers et al., 2015). 


As athermal hydraulics example of the processes that are occurring, consider a salt at an interme- 
diate temperature bounded between a hotter and colder wall (assumed to be ‘gray’). The emission 
spectrum of the hot wall, at temperature Ty, will be in the shorter wavelength region (the black 


12 There are several definitions of the spectral banding of IR. ISO 20473:2007, Optics and photonics — Spectral bands defines 
the visible spectrum to be 0.38um to 0.78 um, near-IR to be 0.78um to 3m, and mid-IR to be 3m to 50yum. Other 
disciplines such as astronomy or lighting define the banding differently. 

13 A correlation and coefficients for the UV absorption (the Urbach tail region) for the same alkali halides considered by Li 
(1980) (and more) can be found, if necessary, in Senatore et al. (1984). 
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Figure 2.5: Thermal radiation emission spectrum at reactor relevant temperatures, com- 
pared to the solar spectrum, and the IR absorption spectrum of water (Palmer and 


Williams, 1974 and Wieliczka et al., 1989). 
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Figure 2.6: Representative absorption spectrum of an alkali halide crystal. Adapted from Li (1980). 
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Figure 2.7: Representative IR absorption spectra of some salts compared to water. The 
500°C and 1000°C emission spectra from Figure 2.5 are included for comparison. 

Data sources: Water: Palmer and Williams (1974) and Wieliczka et al. (1989); FLiBe: Derdeyn et al. 
(2018); FLiNaK: Chaleff et a/. (2016); LiF and NaCl (at 600°C, crystalline solid) Li (1980). 


curves being further to the left on Figure 2.7) than that of the salt at 75. The salt will absorb the 
different wavelengths emitted from the Tw spectrum to different degrees according to its absorp- 
tion spectrum as the IR passes through it, and it will also emit according to the combination of this 
absorption spectrum (recalling that a, = €,) and the blackbody spectrum at Ts (k)Ip)). The cold 
wall will receive a combination of the radiation that has been emitted from the hot wall, plus that 
which has been emitted from the salt, minus that from both sources which has been absorbed in 
the fluid. 


In general the IR absorption spectrum of a salt itself, x,, can vary with temperature, but not sub- 
stantially. It is not even thought to be significantly different between its solid and molten phases: 


« Li (1980) notes that, except for LiF, where a change in absorption was noted when molten, 
compared to when solid, there was no notable variation in molten samples compared to 
crystals (although the molten data was obtained for temperatures that were not substantially 
higher than the melting points). 


Tetreault-Friend et al. (2017) notes that the general behaviour of their molten solar salt ab- 
sorption spectra agree with the Lorentz model, and that the model is considered valid for 
both solids and liquids. 


« There is some dependence of absorption spectrum on temperature, but this is less in amor- 
phous material (molten salts or non-crystalline solids) than for crystal lattices (Modest, 2013, 
Section 3.5). The IR absorption edge moves towards shorter wavelengths with increasing 
temperature. However, with increasing temperature the wavelength spectrum of radiation 
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also moves to shorter wavelengths. 


The presence of fissile actinides and a range of fission products from across the periodic table 
is possible in a liquid fuelled reactor. There is significant information available in the chemistry 
literature regarding the IR absorption of relevant compounds containing many of these elements 
(Mamiya, 1965 for example). This is because IR absorption and Raman spectroscopy are tools 
used in inorganic chemistry to study the structure and composition of these materials. It would, 
however, require some degree of chemistry or spectroscopy domain knowledge to be able to as- 
semble this literature into data to be used in a participating medium heat transfer prediction. No 
readily applicable methods have been found to robustly approximate and combine the absorption 
data for molten salt mixtures, let alone add arbitrary concentrations of other elements. There is, 
however, sufficient evidence that actinides, lanthanides and transition metals can contribute sig- 
nificantly to absorption in the near-IR and visible region, such that the effect cannot be ignored. 
There is a link to redox control in this area, because IR spectroscopy could be used as a method 
for determining the composition and contaminants in a salt (Zhang et al., 2018). 


The conclusion from this high-level assessment of the emission and absorption spectra is that ther- 
mal radiation is important for surface-to-fluid heat transfer, but, because there also is substantial 
potential transparency, it will also be important for surface-to-surface heat transfer. The importance 
will depend on the dimensions of the heat transfer channels and temperatures involved. It is likely 
to be the case that, in many heat transfer analyses, justifying the omission of the modelling of radi- 
ation may require more effort than including it in some form. Given the requirement for salt specific 
optical properties, the lack of comprehensive data, and their variation with composition and con- 
tamination, judging how much effort to put into the thermal radiation component of an analysis will 
benefit from an assessment of its significance to the figures of merit (via SA and UQ). 


Scattering 


The available literature generally does not consider salts to be scattering media to a great extent, 
at least in the absence of significant impurities and contaminants. Highly contaminated salts may 
need to be considered as scattering media, and methods exist to model this, but the increase in 
absorption that may arise from the necessary level of impurity may dominate any scattering that 
arises. 


Refractive Index 


The refractive index of a medium, n, plays a role in thermal radiation heat transfer in semi-transparent 
media. A number of sources of values for relevant salts are available, for example: 


« From Modest (2013, Figure 3.16), n = 1.3 for NaF and n = 1.55 for NaCl or KBr as room 
temperature solids at near-IR wavelengths. 


+ Li (1980, Table 1) has values for optical dielectric constant €., ~ n* for LiF (n = 1.38), NaF 
(n = 1.08) and NaCl (n = 1.53). 


For comparison, n ~ 1.3 for water at IR wavelengths (Wieliczka et al., 1989), which is similar to the 
example values for salts above. This reinforces the point that the reason that thermal radiation is 
not considered in thermal hydraulic analyses related to water, but is considered in salts, is related 
to the higher absorption coefficient in water at relevant IR wavelengths. The absorption coefficient 
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is related to the imaginary part, k, of the complex refractive index m = n-+ /k, and not to the real 
part, n, which determines the speed of propagation (or equally the wavelength) of the photons. 
It should be noted that n is present in Planck’s law, and so influences the blackbody emission; 
the example spectra plotted at 500°C and 1000°C on Figure 2.5 are shown with n= 1.4 asa 
representative value. 


Neutronics 


Energy in fission reactors is generated by the splitting of atoms of plutonium and uranium. As the 
atoms split, energy and neutrons are released along with fission products. The neutrons released 
promote more fission events to occur, which in-turn generate more neutrons, energy and fission 
products. The resultant power generated by the neutron reactions is directly proportion to the 
quantity of fissile material in the core region and the flux of neutrons. 


While neutronics calculations are necessary for all reactor types, the topic is described in this 
volume because there are some features specific to MSRs, and also there are similarities in the 
concepts and modelling approaches between participating medium radiation and neutronics. 


With a sufficient mass of fuel and the correct positioning of devices, such as control rods or mod- 
erators, the reactor reaches criticality, where there is a balance between the neutrons generated, 
absorbed, scattered and lost. This balance is quantified by the effective neutron multiplication fac- 
tor, generally known as the k effective or ka. The balanced (critical) state is achieved when the 
value of ke is around 1. With values below 1 the reactor is subcritical, while with values well in 
excess of 1, the reactor is supercritical, resulting in an exponential increase in the number of neu- 
trons and loss of control of the reactor if no action is taken, or physical feedback mechanisms do 
not sufficiently reduce multiplication. 


The aim for safe operation of nuclear reactors is to maintain a value of keg as close to 1 as possible, 
to ensure that the nuclear chain reaction is self-sustaining, and the aim of neutronics calculations 
is to predict how to achieve this. From the beginning of the use of nuclear reactors, numerical 
solvers for neutronics were developed with the aim of determining the distribution of the flux of 
neutrons in the reactor core. This allows the prediction of the quantities of fuel required to reach and 
maintain criticality, the rate of fuel use, and the effectiveness of active or passive control measures 
to maintain safe operation or shutdown the reactors. 


Accurate predictions of the neutronic state of a reactor also enables accurate predictions of the 
burn-up of the fuel, which in turn allows operators to account for the material inventories of the 
reactor: 


Two important analyses involved in licensing any reactor design are those of mass ac- 
countancy and radiological source terms. Mass accountancy calculations are typically 
done from a nuclear safeguards and non-proliferation perspective. The goal of these 
calculations is to confirm that from fuel production to disposal, the exact amount and 
location of all nuclear material is both knowable and known. While these calculations 
require accounting for a broad scope of possible factors, one of these factors is the 
evolution of the nuclear material while in the reactor. Radiological source term analy- 
sis, on the other hand, focuses on the safety implications of a particular design. This 
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analysis asks what the total radiological source is in terms of possible dose rates to 
the surrounding areas in the event of an accident scenario. 


Graham et al., 2021 
These issues are particularly important in liquid fuelled reactors, where the nuclear material is not 
held stationary and confined in solid fuel elements, and are discussed further in Section 4.2. 
Cross Sections and Energies 


The nuclear reactions that are expected in a reactor are quantitatively described in databases of 
neutron ‘cross sections’, which represent the probability of a neutron with specific incident energy 
interacting with a given isotope. 


The neutron cross sections used in reactor calculations are held in data libraries, such as the 
ENDF/B-VII.1, ENDF/B-VIII.0 or JEFF 3.x libraries that are available from Brookhaven National 
Laboratory'* or the NEA’. 


105 - 
= | 
5 | 
2 4,0 
= 10 
2 F 
Oo 
o | 
wn 
S | 
3 L 
Op ac 
10° | 
if 235 total 
- 2851) fission 
i 285) capture 
pee el eel 
10°12 101° 10° 10% 107 10? 10° 10° 
10° Thermal Resonance Fast 
= | 
5 | 
2 ino 
= 10 
2 r 
o 
9 L 
n 
oi L 
wn 
o L 
0 44-5 
10° 
r 298 fission 
r 238 capture 
bee el el 
fo? 197° 10° 10° 107 10? 10° 10° 


Energy (MeV) 
Figure 2.8: Comparisons of the microscopic neutron cross sections of the fissile (7°°U) 


and fertile (°U) isotopes of uranium at 21°C. NEA JEFF-3.3 library data. The cross 
section unit is 1 barn = 10778 m?. 


The microscopic neutron cross sections for the fissile (7°°U) and fertile (7°°U) isotopes of uranium 
used in nuclear fuels are plotted in Figure 2.8 as an example. The reactions with available cross 


14 Brown et al. (2018), www.nndc.bnI.gov/exfor/endf00.jsp 
15 Plompen et al., 2020, www.oecd-nea.org/dbdata/jeff/jeff33/index.html 
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section data are the fission of nuclei, the capture of neutrons by nuclei and the scattering (elastic 
and inelastic, contributing to moderation) of neutrons by nuclei they encounter. Three categories of 
neutron energy are indicated on the figure: thermal neutrons, with energies of approximately 1 eV 
or less'®, resonance neutrons associated with the region of many peaks and fast neutrons with 
energies of 1 MeV to 20 MeV. 


The microscopic cross section, denoted by o, characterises the probability of an interaction of a 
neutron with an individual nucleus, and has the units of area. Multiplying this by NV, the number 
density of a given type of atom (units 1/m%), gives the macroscopic cross section £ = No (units 
1/m). This represents the probability per unit path length that neutrons will interact with a nucleus 
in a region of material of a given composition. 


Figure 2.8 also demonstrates the different fission and capture behaviour between *°°U and 7°°U. 
Fission reactions dominate over capture reactions for 7°°U, while 7°°U has the opposite behaviour. 


7 of the nuclide in- 


A large range of factors influence the cross sections. For example, the ‘spin 
volved in the scattering interactions plays a significant role in the degree of scattering the neutrons 
undergo. The collision trajectory of a neutron’s motion relative to the nuclei will influence the prob- 
ability of whether it ‘bounces off’ a nucleus and scatters, or is captured. The nuclide temperature 
influences the cross sections through ‘Doppler broadening, where the resonance absorption peaks 
broaden due to the increase in the thermal motion of nuclei at higher temperatures. An introduction 


to these effects can be found in Duderstadt and Hamilton (1976). 


For calculations, the neutron cross sections can be applied either in a continuous way (every data- 
point is considered for each nuclide) or they are divided into a number of groups of differing ener- 
gies. The distribution of the groups depends on where on the energy spectrum of incident neutrons 
the reactor is considered to operate. For example, LWRs are thermal spectrum reactors whose 
neutron cross section can be adequately approximated by two energy groups, while fast spectrum 
reactors will require more energy groups to approximate the distribution of absorbance peaks in 
the resonance region of the spectrum. The discretisation of the neutron cross sections considered 
for fast reactors is dependent on the numerical method used and the geometry considered. 


Neutron Transport 


The different equations that are solved by the variety of methods applied in modelling the neutronics 
are derived from the Boltzmann equation for the transport of neutrons (also called the neutron 
transport equation). 


The neutron density, depending on the location in space, angle (direction) of travel, energy'® of the 
neutron, and varying in time is denoted by 


n(r, Q, E, t) 


This means that, at time t, the number of neutrons in an infinitesimal volume d?r around the 


16 9.025 eV at room temperature, corresponding to a mean speed of approximately 2200 m/s 

17 Here spin refers to the intrinsic angular momentum carried by elementary particles (quarks, leptons and bosons), hadrons 
(neutrons and protons) and atomic nuclei. Neutrons and protons have a spin of 1/2, so the addition of a neutron, forming 
a different isotope, could change the spin state of a nucleus and therefore it could affect the scattering behaviour of the 
isotope. 


18 In the form of kinetic energy, hence related to speed, E = $™n v2 
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position r, travelling a direction "9 contained in the infinitesimal solid angle d?2Q and an energy, 
E, contained in the infinitesimal energy dE is 


n(r, Q, E, t)d?rd?QdE 


The angular neutron flux, used in the neutron transport equation, is obtained by multiplying the 
angular neutron density by velocity 


and integrating this over all possible angles of neutron travel gives the scalar neutron flux 
o(r, E,t) = [uw Q, E, t)d?Q 
An 
The flux of neutrons of all energies is given by 
g(r, t) = [oe E,t)dE 
0 


Terms using ~ and ¢ are combined to give the neutron transport equation. The description of the 
quantities and the terms above and the equation introduced below are a summary of the step- 
by-step derivation that is presented by Duderstadt and Hamilton (1976) or Demaziére (2020), and 
it is recommended that these references are consulted for more detail, and further guidance on 
approaches to numerical solution. 


1 0 x ss r . 
V(E) Ot (r,Q,£,t) + Q-Vyp(r,Q,E, t) + X2(r, E, thp(r, Q, E, t) 


=f fr. E' > £,Q! = Q, t)b(r, EQ, t)e?Q' dE’ 
4n 0 
E) 


co i | 
ee fo 7 Br, EN) (ENE e(r, Ee, t)d(r, ag t)dE" mu = xa) ci(r, t) 
0 i= 


The macroscopic neutron cross sections are introduced to the neutron transport equation by X;, 
x, and *&,, for all possible collision reactions/interactions (total, t), all fission interactions (f) and 
double-differential@! scattering (s). 


The average number of neutrons produced per fission is represented by v. A fraction of (1 — 6) of 
these are released as ‘prompt’ neutrons by a fissioning nucleus, and the rest are ‘delayed’ neutrons, 
originating from fission products that have a relatively short half-life (seconds to a minute), and emit 
a neutron when they decay. 


The probability density function for the energy of neutrons resulting from fission events is x,. The 


19 Where {2 is the direction unit vector, usually expressed in spherical coordintates, related to the neutron velocity vector by 
Q =v/|vI. 

20 Meaning simply that it depends on angle. 

21 ‘double’ in the sense that it includes angle and energy changes simultaneously. 
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delayed neutron precursors are accounted for in Ng groups, and the neutrons introduced by the i” 
group are emitted with probability density function xq,; for their the energy. The DNPs are grouped 
by how long they take to decay, with a corresponding decay constant »;. The concentration or 
number density of DNPs in each group is C;. There are usually six (Sometimes eight) DNP groups 
considered. The fraction of prompt neutrons is related to the fractions in the delayed groups by 


Na 
BE) => Blt) 
i=1 


The neutron transport equation is stated above in integro-differential form, and despite looking 
complex, is only counting the neutrons entering and leaving an arbitrary volume, V, enclosed by a 
surface, S, that are characterised by a specific energy, F, and are travelling in a specific direction, 
Q. This equation has the terms 


D+ @+@=-O©+O+ © 


a) The rate change of neutrons contained in V. 
(2) The net change of neutrons streaming across S. 


(3) The loss of neutrons at this specific E and Q from V as a result of all collision reactions. For 
example, a neutron scattering and only changing direction and/or energy moves to a different 
EorQ. 


(4) This is the inscattering term which represents neutrons that come from a different E’ and 92’ 
scattering into this E and Q. This term is an integral because it accounts for the neutrons 
scattering in from all other angles and energies. 


(5) The emission of prompt neutrons, depending on the scalar flux of neutrons at a different 
energy E’ inducing fissions that release neutrons in the energy in question E. 


(6) The emission of delayed neutrons from each group. 


The DNP concentrations C; need to be calculated using Ng transport equations of the form 


co 


a t) +V- UC; (r, t) —V.- r;VCi(r, t) — fic E)v(E)X-¢(r, E, t)¢(r, E, t)dE - NiCi(r, t) 
ee ———_ —___ 
convective and diffusive transport 0 


where the first term on the right hand side is similar to (5), in that it is the rate of production of 
precursors by fission that will lead to the delayed emission of neutrons. The second term is the 
rate of precursor decay. The need for the convective and diffusive transport terms is specific to 
MSRs because the DNPs are mobile, while most neutronics calculations are based on solid fuels 
held in either fuel assemblies or in TRISO particles. U is local flow velocity vector of the salt and I; 
is the appropriate diffusivity of C;, taking molecular diffusion and turbulent transport into account 
(although this information may not be readily available in simplified flow models). 


This mobility of the fuel means that the neutronics is strongly dependent on the local fluid density, 
with reduced reactivity in regions of higher temperature (hence correspondingly lower density). 
Doppler broadening also depends on local temperature, and results in a decrease in the reactivity 
for most reactors, although this usually has a smaller influence than density changes. Despite 
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their contribution to overall neutron balance being relatively small, the motion of DNPs can lead to 
problems with reactor control, such as the way that flow rate changes can cause reactivity changes 
for the reasons discussed in Section 2.1. 


The neutron transport equation has 7 independent variables (3 for r, E, 2 for Q, t) plus Ng aux- 
iliary C; equations, and while it is the conceptually most simple description of the processes, it 
cannot be readily solved. It requires approximations to be made and numerical methods applied; 
these have sophisticated underpinnings and steps, which can be complex to implement. The need 
for convective DNP transport also requires the flow solution to be coupled to the neutronics. Sec- 
tion 3.4 discusses the various calculation methods for neutronics and the coupling methods that 
are available. 
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This section describes the methods and tools that are available to model single-phase flow and 
heat transfer in geometries and flow conditions relevant to molten salt reactors, as well as giving an 
overview of the methods available for modelling neutronics. This scope does not cover the details 
of modelling dissolved gas retention and transport, species transport and inventory tracking, or 
solidification, although some description of how these topics can be coupled to thermal hydraulics 
is discussed in Section 4.2. 


Convection Heat Transfer Correlations 


Heat transfer by convection in molten salts is not fundamentally different to that in other fluids (as 
discussed in the reviews described in Section 3.1.1 below). Salts have a higher Prandtl number 
than gases or water, and lower Reynolds number turbulent flows are expected, but neither of these 
features is intrinsically difficult to model. Internal flow heat transfer correlations that are valid for Pr 
values relevant to molten salts are well established. The recommended heat transfer correlation for 
fully developed turbulent flow', applicable to constant wall heat flux and constant wall temperature 
conditions, is the Gnielinski correlation: 


(f /8)(Re — 1000) Pr 
Nu= 
1 + 12.7(F/8)1/2( Pr?/5 — 1) 


which has a validity of 2300 < Re < 5 x 10° and 0.5 < Pr < 2000. This is a correlation for a 
smooth circular pipe, where the friction factor is evaluated using 


4 
(1.58 In(Re) — 3.28)? 


It also possible to include the effect of surface roughness, which tends to increase heat transfer, 
directly through f by using a correlation such as that of Haaland: 


4 


3.4735 — 1.5635In | ( 2 ere 63.635 \ | 
: : D Re 


which has a validity of 4000 < Re < 10° and 2e/D < 0.1. This approach is only an approximation 


for the effect on Nu (Incropera et al., 2011). Correlations specifically designed for fully rough? 


This means both hydrodynamically and thermally fully developed. References such as Kakag et al. (1987) and Rohsenow 
et al. (1998) contain the correlations for Nu and f shown here, and also describe methods and correlations for assessing 
Nusselt numbers in specific developing flows, such as entrance regions. 

Roughness Reynolds number, Ree = eu,/v, using the roughness height ¢, is used to characterise a surface as being 
either hydrodynamically smooth, Re. < 5, transitional, 5 << Ree < 70, or fully rough Ree > 70. 
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channels are documented in Kakag ef a/. (1987, Chapter 4), which also provides the means to 
modify any smooth pipe Nusselt number correlation, Nus, for the effects of transition roughness by 
using the ratio of rough to smooth friction factor: 


fe (£)’ sare n=0.68Pr°*° if1< Pr<6 
Nu, \f neil if Pr > 6 


This holds for f/f; < 4; for larger friction factor ratios than approximately 4, there is no further 
increase in Nu. The relationship for n indicates that convection for fluids with Pr > 1 is more 
affected by surface roughness than, for example, gases, where Pr ~ 0.7. 


The frequently applied Dittus-Boelter equation is not recommended for flows where Re < 10* and 
is also less accurate than the Gnielinski correlation at higher Reynolds numbers (Kakag et al., 
1987, Chapter 4). Tube bundle heat transfer correlations for cross flow and axial flow have similar 
validity ranges for Pr > 1000, and are also not inherently limited by the properties of salt (Kakag 
et al., 1987, Chapters 6 and 7). 


Lower Re flows and a reliance on natural circulation in some situations does mean that there will be 
a requirement to assess whether a correlation suitable for natural or mixed convection is needed. 
Free convection is well documented, but mixed convection less so (see Volume 3). Jackson et al. 
(1989) provide a thorough review from an experimental and theoretical point of view for ‘moderate 
Pr fluids’ (not liquid metals) in vertical tubes. It is worthwhile restating the review’s summary: 


In laminar mixed convection, heat transfer is enhanced in heated upward flow and 
impaired in heated downward flow. The problem is amenable to calculation, although 
transition to turbulent flow may occur earlier than for forced convection or free convec- 
tion alone. 

Turbulent mixed convection, heat transfer to moderate Prandtl number fluids is dictated 
by changes in turbulent diffusion. In heated upward flow heat transfer is impaired with 
modest buoyancy and enhanced with high buoyancy. It is not possible to correlate with 
precision heat transfer in the region of impairment using local parameters. In contrast 
heat transfer levels in heated downflow increase monotonically with increasing buoy- 
ancy and have been correlated successfully in terms of local variables. 


The expected magnitude of these effects should be assessed wherever buoyancy is expected, and 
it should be specifically noted that heated upwards flow? can be complex and hard to predict with 
simple methods in some conditions. 


Although there are not significant conceptual differences to other fluids, and the same modelling 
and correlations can be used in most circumstances as a first approximation, there are nuances 
and differences in emphasis that are important to consider for convection in MSR applications. 


Property uncertainty: while the same correlations may be used as other fluids, their accuracy is 
tied to that of the fluid properties used to evaluate them. k, u, cp and p¢ are all required to 
determine a Heat Transfer Coefficient (HTC), and so uncertainty in properties (either caused 
by lack of data or changes in composition) will propagate to the heat transfer prediction. The 


3 otherwise called ‘buoyancy assisted flow’, which also occurs in downwards flows that are cooled. 
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methods that can be used to assess the effect and significance of uncertainties in material 
properties when applied in a system code or in a CFD calculation are discussed in Volume 4. 


Property variation: the properties of salts depend on temperature, and the effect of this on con- 
vection will be discussed in Section 3.1.2. 


Heat generating fluids: the presence of heat sources in fuel salts alters their convection be- 
haviour, and is discussed in Section 3.1.3 


Simultaneous radiative heat transfer: where salts are semi-transparent there is heat transfer 
between fluid and walls simultaneously with convection. This will be discussed in Section 3.1.4. 


Uncertainties and complexities in the use of correlations that arise from these effects will also need 
to be compared to the inherent uncertainty in any correlation and the situation to which it is applied. 
For example, the assumptions of fully developed flows and constant wall heat flux (or temperature) 
are rarely truly representative in realistic situations. 


Experimental Studies and Correlation Reviews 


Accompanying the resurgence of interest in MSRs, a number of reviews have been conducted 
into the available experimental data and assessments made of the applicability of heat transfer 
correlations: 


Ambrosek et al. (2009) re-evaluate previous convection experiments, including those from ORNL 
related to the MSRE, in the light of improved knowledge of salt material properties. This 
supports the conclusion that convection behaviour in salts is not intrinsically different to other 
fluids, however, the high temperatures involved mean that attention needs to be paid to heat 
transfer fundamentals and scaling analysis in designing and interpreting experiments. The 
combination of convection and radiation will need to be understood, and the need for high- 
quality property data is essential, preferably with an estimate of its uncertainty. Britsch et al. 
(2019) also includes a review of historical experiments that supports these conclusions in its 
description of creating a contemporary convection loop experiment. 


Holcomb and Cetiner (2010) and Yoder (2014) include a summary of historical heat transfer ex- 
periments from ORNL and the applicability of conventional correlations for internal forced 
convection is assessed. This demonstrates that there is not much heat transfer data avail- 
able, and that the difficulty of performing the experiments, as well as the age of some of the 
data (mid 1950s) makes the potential uncertainties relatively large. Conventional heat trans- 
fer correlations are shown to be applicable for laminar and turbulent flow, as long as suitable 
material properties are used. The potential for the formation of high thermal resistance lay- 
ers is also noted in conditions where there is inadequate corrosion control, which can distort 
results. The form of the correlations used is discussed in Section 3.1.2 below. 


Wu et al. (2012) report experiments for turbulent and transition Reynolds numbers for HITEC and 
LiINO3, which are compared to data from the same historical experiments as above, and to 
similar convection correlations that also account for temperature varying properties. Similar 
conclusions to that of Yoder (2014) are reached, where the correlations compare relatively 
well to the data, noting that it contains scatter and a typical agreement of approximately 

+20 % or better is seen. 


Yoder et al. (2018) report experimental results for natural convection from a vertical cylinder in 
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a larger cylindrical space and demonstrate that the natural convection correlations for this 
configuration developed for other fluids can be successfully applied. The flows are predomi- 
nantly laminar, and there are more uncertainties for turbulent results (generated by high Ra 
conditions). It is also noted that there is a need to correct the assessment of convective heat 
transfer for simultaneous radiation. 


Zhang et al. (2020) compare correlations of the kind discussed in this volume, the same experi- 
mental data as discussed above (with attention paid to modifying historical results with im- 
proved knowledge of material properties) and the results of CFD models using coolant and 
fuel salts. The authors are able to find an overall good agreement between the three ap- 


proaches, but note that uncertainties of up to +20 % are seen. 


Temperature Variation in Properties 


Heat transfer correlations for internal flows are typically derived with the assumption of constant 
fluid properties over the cross-section of a channel at a particular location, and the properties are 
usually evaluated at the bulk temperature (as they are in the Gnielinski correlation above). Bulk tem- 
perature is the mass flux weighted average temperature of the fluid, see Volume 2, Section 2.1.2. 
However, this assumption of property uniformity does not hold when the material properties depend 
strongly on temperature, particularly viscosity. In cases where there is significant heat transfer at 
a wall, the viscosity in the near-wall region can be substantially different, leading to changing ve- 
locity gradients, and a change in convection behaviour. This is the case for salts, where, because 
of their relatively high Pr, there can be significant temperature differences across the momentum 
boundary layer. 


Kakag et a/. (1987, Chapter 18) provides the most detailed description of the effects of temperature 
on convection correlations (although Kays et al. (2005) and Incropera et al. (2011) also have some 
useful guidance), and suggests that for liquids* the following corrections can be applied 


Nu _ ( be 
wes = (2) 
f _ [{ Bb 7 
aC) 


The ratio of viscosity evaluated at the bulk (b) and wall (w) temperatures is used to correct the 


Nusselt number evaluated assuming constant properties (CP). Various values to choose for n and 
m are reviewed in Kakag et al. (1987), and recommended values for turbulent flow that can be 
applied to the Gnielinski correlation are 


Me 


Bw 
He 


bow 


>1, n=0.11 for heating liquids 
<1, n=0.25 for cooling liquids 
which is valid for 10* < Re < 5 x 10®, 2 < Pr < 140 (Re and Pr evaluated at the bulk temperature) 


and 0.08 < tw/, < 40. As an example, for FLiBe at 850K bulk and 950K wall temperatures, 
the viscosity ratio is 1.59 and ratio of ratio Nu/Nucp = 1.053 using n = 0.11, meaning that heat 


The behaviour for gases is different — viscosity increases with temperature, and density and conductivity also change 
significantly. For liquids, viscosity tends to decrease with increasing temperature, and varies more than other properties. 
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transfer is expected to be increased by approximately 5% due to the variation of properties. 


An alternative correlation accounting for property variation in smooth pipes, that does not use the 
viscosity ratio approach, and is valid for liquids and gases is that of Sleicher and Rouse (1975) 


Nup = 5 + 0.015Re?Pr?, 
a = 0.88 — 0.24/(4 + Prw) 
b= 1/3 + 0.5e-0-6Prw 


valid for 10* < Re < 10° and 0.1 < Pr < 10°. It also a good example of the need for care and 
clarity in which temperature is used to evaluate properties: Nu,, Res and Pr, use the bulk, film 
and wall temperatures respectively in the same correlation. It is notable that turbulent correlations 
with validity between transition Reynolds numbers and Re = 10* are unusual — the Gnielinski 
correlation was created to specifically extend to lower values. 


For laminar flows the Sieder-Tate correlation can be used, which includes a term where n = 0.14 


1/3 0.14 
Nu = 1.86 (Rerr? ) (=) 
L Hw 


and there are some features of this correlation that raise points relevant to the reviews mentioned 
in Section 3.1.1 above: 


* This correlation contains the Graetz number term, Gz = RePrD/L, making it dependent 
on downstream distance, so is necessarily only applicable to flows that are still developing 
(which in laminar conditions can be over a long distance). 


« The correlation is derived for constant temperature boundary conditions — for laminar flows, 
the difference between constant temperature and constant heat flux boundaries matters; for 
turbulent flows, often the same correlation can be used by both. 


« The value of Nu that this correlates is the mean Nusselt number over the entire channel 
length, L —it does not give the local value varying as a function of L/D. This also introduces 
some dependency on the upstream and inlet conditions to the channel, and the correlation 
has limits on its applicability, making it not accurate for very short and very long channels. 


* The values of Nu, Re and Pr relate to material properties at the bulk temperature. 


When the correlations are reviewed in the contemporary references it is often the case that the 
descriptions are incomplete with regard to these points, and the guidance is somewhat ambigu- 
ous. The original references or thorough assessments such as in Kakag et al. (1987) should be 
consulted. 


The experiment and analysis by Cooke and Cox (1973) of forced convection heat transfer using a 
proposed MSBR fuel salt, Lip-BeF,—ThF,—UF,, is a good example. They show that the laminar 
Sieder-Tate correlation above, and its fully developed turbulent equivalent (also applying n = 0.14) 
agree relatively well with the data when correctly compared (Figure 3.1). They report that the corre- 
lations lie about 13 % above the experimental data in the turbulent region, and also propose slightly 
modified coefficients to the relations that they use that bring the agreement to within approximately 
6%. In this reference, the need for accurate properties is emphasised, which is true in all instances, 
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Figure 3.1: Forced convection data compared to correlations accounting for property 
variation across the laminar, transisiton and turbulent flow regime, from Cooke and Cox 
(1973). ws is the viscosity at the surface temperature and y is at the bulk temperature. 


but this also shows that salts do not behave in a different way to other fluids for forced convection 
heat transfer. 


The concluding remarks in Kakag et a/. (1987, Chapter 18) discuss how more detailed analytical 
assessment is not tractable, requiring numerical methods, and that experiments are hard to perform 
because large wall-to-fluid temperature differences are needed. It is also the case that (minor) 
improvements in accuracy are possible for turbulent flows by considering the specific shape of a 
non-circular channel (rather than a simple circular channel equivalent hydraulic diameter). 


This suggests that CFD assessment of the geometry, flow and heat transfer regime and materials 
specific to a problem at hand would be a valuable way of gaining insight and developing correla- 
tions. This requires that the chosen CFD turbulence modelling approach is able to represent the 
interaction between property changes and turbulence (as well as the effects of buoyancy and heat 
generation, as discussed below). 


An example of determining the effective heat transfer and friction factor from CFD is given in 
Study B: Fuel Assembly CFD and UQ for a Molten Salt Reactor. 
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Natural Convection 


Natural convection (free convection and buoyancy driven flow in enclosures) inherently takes den- 
sity differences into account, although sometimes only as the Boussinesq approximation, but in 
many cases correlations do not account for other property changes. Kakag et al. (1987, Chapter 
18) again provides a detailed assessment, and Yoder et a/. (2018) reviews similar correlations. 
Corrections that involve the ratio 4../Hw are applied and the need to incorporate a temperature 
varying value of 6 is noted, where it appears in the Grashof and Rayleigh numbers. The experi- 
ments and numerical modelling of Cuckovic-Dzodzo et al. (1999) for laminar recirculation in cavities 
for high Pr fluids shows that accounting for property variation is necessary to correctly represent 
the flow patterns and heat transfer. 


In some cases, improved natural convection correlations are achieved by modifying the tempera- 
ture used to evaluate properties that parameterise them. Reference temperatures of the form 


T, = Tw — X(Tw + Too) 


with values of 0.25 < x < 0.5 may give the best result; x = 0.5 gives the familiar film temperature 
that is typically used for natural convection. 


One of the conclusions that can be reached by reviewing the limited experimental data, and the 
careful analysis that accompanied it in historical references, is that it may serve its most useful 
purpose as guidance on interpreting, parameterising and generalising high quality CFD results 
into correlations, rather than as a source of correlations itself. 


Fluids With Internal Heat Generation 


Heat generation within a fluid means that the temperature distribution implicit in the heat transfer 
correlations, derived from idealised constant wall heat flux or constant wall temperature boundary 
conditions, is not applicable. This is particularly the case for laminar flow, where the shape of 
the temperature profile in a channel has a substantial effect. Typically turbulent correlations are 
less sensitive to the details of a channel cross-section geometry (and so they can usually be 
successfully parameterised by a hydraulic diameter). However, if the volumetric heat source (qyoy, 
in W/m?) is sufficient, it will change the temperature distribution in turbulent flows too, increasing 
temperatures particularly in the slow moving region near to a wall. For a given wall heat flux into a 
channel from the wall (heating the fluid), the presence of heat generation in the fluid will increase 
the wall-to-bulk temperature difference, and reduce the apparent HTC. For wall heat flux out of 
a channel (cooling the fluid), the opposite is true, wall-to-bulk temperature difference is reduced, 
indicating a higher apparent HTC. 


Heat transfer from flows with internal generation was successfully analysed for hydrodynamically 
fully developed® flows in the pre-MSRE era (Poppendiek and Palmer, 1952, Poppendiek and 
Palmer, 1955 and Siegel and Sparrow, 1959). Subsequent experimental studies showed that the 
predictions made by the analysis agreed well with experiments for simple pipe flow configurations. 
For example, Kinney and Sparrow (1966) show experimental wall-to-bulk temperature difference 


i.e. the velocity field was developed in a long adiabatic section with no heat source, but the temperature field does not need 
to be fully developed. 
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results for turbulent flows of a fluid with Pr = 3 to 4 that agree with the analysis. In this case, 
it was demonstrated that, as predicted, when a pipe is adiabatic and there is a volumetric heat 
source in the fluid, the wall temperature is actually higher than the bulk (or centreline) temperature. 
Michiyoshi et al. (1968) shows experimental data for temperature profiles across a channel for 
turbulent flows, which compares well with the analysis of Siegel and Sparrow (1959). For laminar 
flows in vertical pipes, they demonstrate that buoyancy effects near to the wall would need to be 
taken into account in the analysis. The validation was largely intended to test the applicability of 
Reynolds analogy in heat generating fluids — the hypothesis that momentum and thermal diffusiv- 
ity are similar in moderate to high Pr fluids was used as an assumption in the derivation of the 
predictions, and confirmed by the experiments. 


More recent analysis using a similar approach (Di Marcello et a/., 2010) was also able to predict 
the same temperature profiles, and then express a relationship (Fiorina et al., 2014) between 
the turbulent Nusselt number for fluids with internal heat sources, Nug,,,, and without, Nu (i.e. a 
conventional pipe flow correlation). This is a function of heat source, q,o/, wall heat flux, q (signed 
as inwards or outward; inwards is positive), pipe diameter, Re and Pr: 


Nu ger pom, 1 


Nu qvoiD _ 1.656 
ik ar q Pr®-4 Reo-5 


Therefore, there is a large difference in apparent Nusselt number expected for large values of D, 
low values of Pr or Re, or large ratios of qyoi/q. This expression is not necessarily convenient to 
implement in a conventional system code, because it contains q, the pipe wall heat transfer, which 
is usually the quantity that a Nusselt number correlation is being used to evaluate in the code. 
Therefore, there would be a need to iterate to find the suitable value, and a code would need to be 
written to expect this. To implement and test functionality for wall heat transfer with heat generating 
fluids in a calculation, the 1950s references above are recommended as a starting point, because 
of the detail and clarity of their descriptions. 


The analysis that is discussed in the references above for turbulent flow is for axisymmetric, fully 
developed flow in a pipe, assuming constant fluid properties, uniform heat source, and not con- 
sidering buoyancy effects, so the ability to apply it to realistic conditions has some limits. In these 
analytical approaches, the success of the prediction is dependent on knowing, as an input, an ac- 
curate distribution of velocity and eddy diffusivity across the channel. This is available for simple 
configurations like fully developed flow in smooth circular pipes, but is not known in general. How- 
ever, the conclusion that, with an appropriate description of turbulence, the presence of internal 
generation is readily predicted supports the application of CFD to these conditions. CFD using the 
material properties, geometry and heat transfer conditions for a specific situation would therefore 
be more insightful, both as an analysis method in itself, and to generate inputs and correlations for 
system code models. The references discussed here would be valuable for the interpretation of the 
results, and to help define the appropriate parametrisations to extract transferable correlations. 


The situation of a heat generating fluid in a closed vertical cylinder (i.e. no net flow through a 
channel) is studied in detail in Study B. 
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Combined Convection and Radiation 


Thermal radiation exchange between a scattering®, absorbing and emitting flow and a wall has 
been a relatively well studied topic from the 1960s (Einstein, 1963) onwards, because it has ap- 
plications in combustion, high-speed aerodynamics (where hypersonic shockwaves can radiate 
significantly) and in nuclear rocket propulsion. 


The main consideration is a fluid’s absorption coefficient, «, and the consequent optical thickness 
tT = KL of the fluid in a channel of dimension L (which would be the radius for a circular channel). 
For the condition of a hot channel wall heating an internal flow, there is an increase in apparent 
total heat transfer and Nusselt number caused by the participating fluid. The enhancement to heat 
transfer is optimal at intermediate optical thicknesses: 


« When optical thickness is small (t < 1, optically ‘thin’), the fluid does not absorb energy, so 
radiation does not contribute to cooling the wall’. 


* At intermediate optical thicknesses the maximum radiation interaction occurs with absorption 
occurring across the channel width. 


« When optical thickness is large (7 >> 1, optically ‘thick’) then the radiation emitted by the 
surface is absorbed in a thin layer near to the surface and partially re-radiated back towards 
the wall, shielding the central regions of a channel from the direct radiation. 


The mathematical formulation of coupled convection and radiation equations is relatively, complex, 
even for simple geometries, and, like for heat generating fluids, is generally presented for hydrody- 
namically developed flows with constant properties (meaning that the velocity and turbulent dissi- 
pation is needed as an input). The main references for this topic are Kakag et al. (1987, Chapter 
19), Modest (2013, Chapter 22) and Howell et a/. (2016, Chapter 13), and they each present the 
concepts in detail and review the studies in the wider literature. 


Two important non-dimensional parameters emerge in the analysis of combined radiation and con- 
vection: 


The conduction-to-radiation parameter, also known as Stark number: 


kk 


eo 
CR A4noT3 


where a smaller value of Ncr means that situation is more influenced by radiation. Despite having 
‘conduction’ in its name, this parameter occurs in the non-dimensionalised equations for radiation 
combined with both conduction and convection. Lower values have the effect of increasing the size 
of the thermal boundary layer, in a similar way to lower Pr. Convection dominated conditions occur 
with larger values of Ncr. 


The convection-to-radiation parameter, usually called Boltzmann number: 


pUCp 


B — 
. roT3 


In the rest of this discussion, it is assumed that scattering is not significant (as is the case in salts), but the concepts and 
references readily encompass it. 

This is for the case where all of the walls bounding a channel are assumed to be at a uniform temperature, in the case of a 
difference in wall temperatures over a cross-section, what is not absorbed by one section of wall travels across the channel 
and there is radiative heat transfer with other surfaces. 
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Similarly, smaller Bo means more radiation influence, and convection dominates at larger Bo. 


These parameters do not directly give the quantitative ratio of thermal radiation heat flux to that 
from other the heat transfer mechanism, but they are used to qualitatively indicate and compare 
their relative magnitudes. The actual heat flux ratio depends on the absolute temperature, wall to 
bulk temperature difference, and optical thickness, and is non-negligible in channel flows (typically 
increasing convective fluxes by approximately 5 to 20%, depending on the flow conditions). It 
is also noted by Modest (2013) that ‘no thermally fully developed conditions can exist for forced 
convection in a flow combined with appreciable thermal radiation’ because of the non-linearity 
of radiation with temperature. It is also the case that treating channel flows as 1D can lead to 
inaccuracies because axial radiation can be significant, especially where channels are converging 
or diverging, changing the angle of (hence of emission from) the walls relative to the channel 
centreline. 


The analysis described in the references above (and also available in a concise form in Chawla 
and Chan, 1980 and Yener and Ozisik, 1986) is able to predict temperature profiles and heat 
transfer, and was generally carried out by numerical solution even in the earliest references. It 
has not, however, been formulated in a way that can be easily incorporated into a heat transfer 
correlation in a system code. Large simplifications can be achieved in such calculations by ignoring 
scattering, variation of properties with temperature and wavelength (spectral dependence), but it 
can lead to inaccurate or unrealistic results (Howell ef a/., 2016, Chapter 13). Even with these 
simplifications, the effect of radiation on convection depends on geometry, material properties, 
flow rate and temperature, and so creating universal correlations would be challenging. The same 
points also apply to predicting natural convection or boundary layer flows involving a participating 
fluid. 


More recently there have been studies conducted by Chaleff et al. (2016), Derdeyn et al. (2018) 
and Coyle et al. (2020) that are specifically related to MSR salts, and in-particular considering the 
variation of absorption for realistic salts with wavelength. 


Given these points, the best approach for gaining insight is likely to come from combining realistic 
geometries, flows and material properties with radiation modelling that is already available in CFD 
codes, as discussed in Section 3.3. The creation of a convection heat transfer correlation that 
incorporates radiation to a participating fluid using CFD is demonstrated in Study B. 


Flow and Convection with CFD 


The topics discussed in Section 3.1 above indicate that CFD is likely to be able to address each of 
the effects individually. Therefore, there is little reason to believe that, wnen combined, the physical 
effects that are most pronounced in MSRs would be expected to lead to a significant reduction in 
accuracy compared to that which would be expected from CFD analysis generally. The modelling 
of turbulence, in combination with unsteadiness and buoyancy, is still likely to be the most chal- 
lenging aspect. However, the presence of thermal radiation and heat sources in near-wall layers 
should motivate the generation of CFD grade validation data, and resolved turbulence simulations 
(Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS)) to assess the performance 
of Reynolds-Averaged Navier-Stokes (RANS) models. 
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The discussion here builds on the CFD guidance in Volumes 1, 2 and 3, and the contrast with the 
more extensive considerations needed for low Pr fluids as described in Volume 5 is also instructive. 
Therefore, only the specific topics of wall heat transfer for moderate Pr fluids and the modelling of 
pebble beds with CFD are discussed in this volume. 


Study B demonstrates an example of applying CFD to a molten salt which simultaneously exhibits 
internal heat sources, strong property variation with temperature, buoyant turbulent flow, property 
uncertainties and participating thermal radiation. 


Near-Wall Modelling 


When modelling wall heat transfer with molten salts, three high level questions arise that motivate 
the rest of the decisions: 


1. What is the effect of the moderately high Pr of the fluid? 
2. Does the thermal boundary layer need to be resolved? 


3. Can wall functions for heat transfer be used, or does the momentum boundary layer need to 
be resolved? 


The answer to the first question is, to a first approximation, that the Prandtl numbers encountered 
are not extreme, and their effect is relatively localised for turbulent flows. Once the distance from 
a wall is large enough that turbulent diffusivity is dominant, then the low contribution of molecular 
contribution makes little difference. The effect is mainly focussed on the laminar sublayer (y* < 5) 
and buffer layer (5 < yt < 30). The boundary layer profiles of Kader (1981), shown in Figure 2.1, 
more than cover the expected Pr range for salts (Kader compares the profiles to experimental data 
with Pr of 5.5, 14.3 and 95). However, they also show that in the viscous sublayer 


Pe Spry 


where 


7+ _ (T-To) 


and 7, = ae 
T, PCpUr 


T; is the ‘friction temperature’ and uses the friction velocity u-. This means that, even when the 
momentum boundary layer mesh is resolved such that y* ~ 1, then more cells would be needed 
to resolve the thermal boundary layer for Pr > 1. This is the second question, and in general, 
the answer is ‘no’ — where the momentum boundary layer is well resolved, there is no analo- 
gous need in general to create a finer mesh that achieves T+ ~ 1. Because the flow is laminar, 
and slow moving at the wall, the heat transfer occurs by molecular conduction. This approach to 
near-wall solution, with a resolved momentum layer, should be able to locally accommodate heat 
generation in the fluid, participating thermal radiation, and temperature dependent properties with- 
out additional modelling considerations, and should be the preferred option where computational 
resources permit. 


There will be some case-specific reasons to test the dependence of the results on first cell T*, 
however. In situations where with high Pr, high heat flux, and strongly temperature dependent 
properties, it may be necessary to check whether a better resolution of the temperature gradient 
changes the results. Another reason might be to assess the onset and development of salt freezing 
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at a surface — for this application Cartland Glover et a/. (2019) resolved meshes to T* ~ 1. 


This leads to the third question — does the momentum boundary layer always need to be resolved, 
or can wall function approaches still be used for heat transfer? The first thing to note is that wall 
functions are generally developed with the assumption of fully developed flows (or equilibrium 
boundary layers), and with the assumption of constant material properties. The latter of these 
two points may itself be sufficient to make the answer ‘no’, conventional wall functions are not 
sufficiently accurate for modelling molten salts where there is a large change in viscosity because 
of local temperature variation. However, it is not the case that the higher Pr inherently means that 
a wall function approach cannot be used. 


Near to a wall, in the laminar sublayer, higher Pr fluids have a higher thermal resistance associated 
with lower conductivity. Thermal boundary layer expressions have ‘P functions’ that account for this, 
and these are used to create wall function models for heat transfer. A common P function is that 
derived by Jayatilleke (1966) from experimental results of fluids with moderately high Pr fluids and 
wall roughness. Malin (1987) reviews a range of P functions and concludes that, other than for 
very high values of Pr (>> 100), then the Jayatilleke P function is suitable. Such models require a 
turbulent Prandtl number, Pr,, but this has been studied quite extensively, as summarised by Kays 
(1994). The temperature profiles of Kader (1981) are also used as a basis for wall functions. 


Recently Sarié and Basara (2018) developed a wall model for higher Pr fluids taking into account 
the variation of Pr; near to a wall. It demonstrated improved performance over the Jayatilleke P 
function models for Pr ~ 20 where the first cell was in the buffer layer, although it was implemented 
in the context of the k-¢-f RANS model, which is based on the elliptic relaxation concept (see 
Volume 3, Section 2.2.4.3), and is not widely available in CFD codes; the approach is also not 
necessarily adaptable to more common models, such as k-w SST. 


In the context of MSRs, Fiorina (2019) shows that wall functions for fluids with internal generation 
need to be adapted. Heat generation in the fluid causes an underestimation of temperature differ- 
ences to a wall when using thermal wall functions, and a modified implementation based on the 
Jayatilleke P function shows good agreement with the analysis and experimental data discussed in 
Section 3.1.3. Participating thermal radiation with high optical thickness, causing significant energy 
transfer in the near-wall layer would be expected to have a similar effect. 


The discussion above was mainly based on the assumption of smooth pipes. Where roughness 
becomes significant, it disrupts the laminar sublayer. For higher Pr fluids, the effect of this on heat 
transfer is more pronounced, because it allows turbulent diffusivity to reach the wall and removes 
the high resistance layer (and exposes the wall to turbulent shear stress, increasing friction factor). 
This effect of roughness is able to be accommodated in wall model boundary conditions in CFD 
codes — the effect of the disruption of the laminar layer is to shift the ‘law of the wall’ velocity 
profile downwards. Even in fully rough® conditions where there is no laminar sublayer, the effect 
of higher Pr is still present to an extent as the final conduction resistance to the wall (Kays et al., 
2005, Chapter 12). For a fully rough wall, some form of wall modelling will inevitably be present, 
so attempting to use a near-wall resolution suited to smooth surface boundary layer that would 
give yt = 1 will not provide the same kind of improved accuracy. Wall modelling, and especially 


8 See footnote 2 on page 41. 
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the implementation of roughness and heat transfer are details that are CFD code specific, and the 
documentation for the approach taken in the particular code being applied should be reviewed. 


Pebble Beds 


Pebble beds represent an application that is well suited to a porous medium representation within 
a CFD code. They have received significant modelling attention in the past (Stainsby ef a/., 2010, 
IAEA, 2013), and the correlations used in system code representations and validation data can be 
readily applied in developing CFD models. 


It has recently become possible to combine CFD grade experimental data (Nguyen et al., 2018) 
and resolved turbulence (DNS and LES, Yildiz et a/., 2020) to provide insight into the details of 
the flows in pebble beds. This is able to be compared to the more conventional correlations for 
drag (such as the well-known Ergun equation) and pebble-fluid heat transfer to create and validate 
porous models of the whole core of pebble bed reactors (Novak ef a/., 2019). 


The references noted here are a starting point that cover (and provide further references to) all 
of the pertinent aspects of modelling pebble beds with CFD. Since this is a rapidly evolving area 
of capability development, the details are not repeated here. Examples of coupling these meth- 
ods to neutronics solvers for the prediction of nuclear heat sources are also discussed briefly in 
Section 3.4.3. 


Radiative Heat Transfer with CFD 


Modelling radiative heat transfer with CFD is introduced in Volume 2, but it does not discuss par- 
ticipating media. Most CFD codes include participating radiation modelling capabilities and couple 
the energy transfer processes in the fluid°. The main choices of approach are the P,, Discrete 
Ordinates (DO) and Monte Carlo (MC) models, all of which can incorporate scattering, absorption 
and emission from a fluid. Two main considerations will guide the choice and cost of a model to 
use to represent radiation: 


¢ What is the optical thickness of the radiation paths in the model? 
* Can the fluid be approximated as being gray? 


Both of these considerations also relate to the state of knowledge of the material and surface 
properties — a costly model is unlikely to be worthwhile if there is poor knowledge of the inputs. 
However, a lack of knowledge about optical properties of a material should not be used as a 
reason to exclude radiation modelling — at high temperatures it can have a significant effect on 
heat transfer in some geometries, and an analysis that does not assess this contribution is likely to 
be challenged. In some cases it can be simpler and faster to include radiation modelling, compared 
to the effort of justifying omitting it. In other cases an unambiguously conservative approximation, 
whose effect on an operational or safety margin is acceptable, can be made by omitting radiation 
in an assessment. 

The energy coupling from the fluid temperature to the thermal radiation field is included, but as pointed out by Howell ef al. 
(2016), radiation and turbulence are coupled. The local temperature fluctuations caused by turbulence are amplified by the 
T* dependence, which increases the emission compared to using only the mean temperature. This effect is included in 


specialised combustion models, where large temperature variations occur, but not in the general purpose models described 
in this section. 
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In general, non-gray behaviour (spectral dependence) is modelled using a ‘banded’ approach, 
where a representative absorption coefficient is applied across a particular wavelength band, and 
the bands are chosen by the user. Adding more bands adds to the cost of the simulation in the P, 
and DO models, because a set of RTEs are solved for each band. For MC, the computational cost 
is associated with tracking particle (photon) paths, and each band will be tracked independently. 
Section 3.3.1 below provides a review of the available guidance on choosing inputs for a banded 
approach. 


The DO and MC methods are able to model any optical thickness, and should be considered as 
being able to provide a reference/comparison calculation for simpler methods. However, this is only 
the case when they are sufficiently well resolved. For DO, the unit sphere is split into segments, 
and higher angular discretisations lead to more cost (solving more transport equations), but without 
sufficient resolution, significant artefacts and inaccuracy can occur, as described by Murthy and 
Mathur (1998) (for the method as implemented in Fluent). MC is a ray-tracing approach, and using 
more rays adds to computational costs, but increases accuracy. 


The P, method is usually recommended (for example by ANSYS, 2020) for optical thicknesses of 
T > 1, and would be expected to be significantly computationally cheaper than either DO or MC. In 
some geometries, using P, for 7 < 1 may be able to provide sufficiently accurate results, but not in 
others (Modest, 2013, Chapter 16), and this should be checked against a more accurate solution. 
It is likely to be the case that using P, will give a more accurate answer than omitting radiation 
altogether, and it is not expected to add significantly to solution costs. 


When optical thicknesses are high (7 > 3), a diffusion based approach (the Rosseland model) may 
be appropriate, although unless there are significant heat sources in the fluid, the difference that 
will be made to convective heat transfer is likely to be small in nearly opaque conditions. Diffusion 
approaches are usually limited to the gray assumption. 


The MC, DO, P; and diffusion approximation methods that are applied to radiative heat transfer are 
conceptually the same as those used for neutron transport, discussed in Section 3.4. P, method is 
the Py method, with N = 1 (related to the number of terms in the spherical harmonic expansion 
that are used); DO is also known as the Sj, method, where AN is related to the number of discrete 
angular directions considered. 


Determining Banded Gray Absorption Coefficients 


Extensive literature is not available describing rigorously transforming the absorption spectrum 
for a semi-transparent fluid, plus additional absorption from specific constituents or contaminants, 
into banded mean absorption coefficients, suitable for inclusion ina CFD model. There is a lack of 
guidance in both how to choose the number of bands, their wavelength ranges, and how to average 
the value of absorption. 


Most of the available literature (such as Modest, 2013) on modelling spectral (non-gray) participat- 
ing medium radiation transport concentrates on the propagation of radiation in molecular gases, 
which has an application in combustion, spectroscopy and atmospheric propagation. For exam- 
ple, the HITRAN database (Gordon et al., 2017) is very detailed, resolving each absorption line — 
millions to hundreds of millions of lines for common molecular gases such as H,O. This data can 
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be assimilated into coarser banded representations, or treated using other sophisticated spectral 
methods. 


The gray approximation is commonly applied in commercial CFD codes, and none of the mod- 
ern spectral methods are available'? — only the weighted-sum-of-gray-gases method, plus some 
models for the effects of particles or soot. For molten salts which have an effectively continuous 
(rather than multiple discrete banded) absorption spectra in the wavelength interval of interest, and 
limited experimental knowledge, these sophisticated spectral approaches are not a ‘good fit’, so 
a simpler approach is required. It should also be noted that, because their intended application is 
radiation transfer in gases, some expressions for participating media omit refractive index, with the 
(sometimes implicit) assumption that it is approximately 1. This is not the case for molten salts. 


Much of the radiative heat transfer analysis for participating media is aimed at predicting the heat 
transfer between an absorbing, emitting and scattering gas and a bounding surface, not for making 
local radiation propagation assessments. It is in this context that the concept of a ‘mean beam 
length’, L,, is introduced, which allows the emissivity of an isothermal participating medium of a 
particular geometry to be related to the absorption coefficient using 


e=1—e “tm, 

Fiveland and Jamaluddin (1991) rearranges this 
_ In(1 —e) 
K=— Ls ‘ 


to produce absorption coefficients for the spectral banded DO model that is available in common 
CFD codes. However, the reference does not provide guidance on the source of the emissivity 
data, or what value of L,, to apply. Similarly an example of this method being applied to produce 
absorption spectra for molten salt (for a solar application) is shown in Tetreault-Friend et a/. (2017), 
where the use of a fixed length is appropriate because of the simplicity of the geometry (a uniform 
layer of material). Furthermore, in ANSYS Fluent, the weighted-sum-of-gray-gases model requires 
a user specified value of L,, or, it will use a default value 
V 
—— 3.67. 
based on the total volume, V, and surface area, A, of the model domain. While these approaches 
may be suitable for uniform layers of material, or for a furnace or combustor (which is typically 
a relatively simple, open domain, with a hot flame region and cooler walls) this is not suitable 
for general purpose nuclear reactor components, where the range of length scales varies from 
tube-spacing in heat exchangers to a large portion of the reactor diameter in an open plenum. 
Therefore, it is necessary to have local, composition (and temperature'') dependent properties 
that are independent of an implicit or assumed geometry, and where an absorption spectrum can 
be split into arbitrary wavelength bands, each of which can have very low or very high absorption, or 
anywhere between. Adapting Equation 11.60 from Modest (2013, Chapter 11.8), the local average 
10 From Chapter 3.3 of Chhabra (2017), written by Michael Modest: ...as of 2015 none of the better known CFD packages 
includes a credible non-gray spectral model. 


11 This flexibility is available in Fluent, for example, where the absorption coefficients in each band can be given a temperature 
dependence via a User Defined Function. 
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absorption is given simply as 


In Modest (2013, Chapter 20.4) the stepwise-gray (or box) model also suggests that this is the 
appropriate approach. 


This is an area of evolving research, and other methods have been suggested. Modest (2013, 
Chapter 11.11) describes the Planck mean absorption coefficient as being a useful quantity in 


evaluating the RTE 
_ i IpykK dr 


EE TO Tay. 


Noting that this is evaluated across the whole wavelength range. (Howell ef a/., 2016, Chapter 
12) describes that applying this as a single value across the whole spectrum can be inaccurate 
in some circumstances, and that other whole spectrum mean definitions are available, but that a 
banded approach offers improvements (Patch, 1967). Coyle et a/. (2020), following the PhD thesis 
of Chaleff (2016), apply what is described as the ‘Planck mean’, but on a wavelength band basis 


_ SR? Inykeadd 

= SR? Tox 

This is not an equation that is evident in, for example, Modest (2013) or Howell ef a/. (2016), and it 
is notable that this approach is not taken in the similar journal paper by Chaleff et a/. (2016), where 
a simple band midpoint is used. It is also not clear whether it is consistent with the definition of the 
banded gray model that it is applied to in Star-CCM+. In the spectral banded DO method defined 
by Fiveland and Jamaluddin (1991), and also stated in the Fluent Theory Guide (ANSYS, 2020) 
and in Modest (2013, Chapter 1.4), the blackbody radiative intensity emission over a band is 


r2 roT* Pe Ep dr 
Ind = (F(nd\2T) — F(n\yT)) ——,  F(ndT) = 2, 
[tora = (FORT) — FOMTIVTLTE, FUT) = Bee 


noting the limits of 0 and oo in the denominator. It is therefore possible that this banded ‘Planck 
mean’ approach is applying double weighting of the spectral variation, and more justification would 
be needed to suggest its use. There are methods discussed by Patch (1967) related to banded 
approximations that make use of ‘group [band] Planck mean absorption coefficients’. They have 
some similarity, but there is significant additional complexity in their subsequent application that 
suggests that their definition is specifically incorporated within those methods. The complexity 
also makes it challenging to readily compare it for consistency with the CFD spectral banded DO 
implementation of Fiveland and Jamaluddin (1991). 


A simple local average is currently recommended. 
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Neutronics 


Existing neutronics solvers are mature software tools, and an overview of the common methods 
and codes is described below. They have assumptions, neutronic data, numerical models and 
implementations based around the heterogeneity of solid fuels. These are likely to perform well 
with reactors that have similar layouts to the MSRE with graphite moderation, but an evaluation of 
the suitability of the neutronic data and the modelling assumptions in them may be necessary for 
fast spectrum MSRs without core internals, such as the MSFR. Similarly, for simulations of LWRs, 
neutronic coupling to thermal hydraulic models is relatively loose because of their static fuel, but 
liquid fuel density changes and the motion of DNP make tighter coupling necessary. 


Therefore, a number of implementations of varying complexity have considered the impact that 
the motion of fluids and the influence of temperature have on the reactor physics of molten salt 
reactors. These measures vary from the modelling of the neutron flux based on one-dimensional 
representation of the fluid motion, to the coupling of fluid dynamics and neutronics solvers (using 
similar coupling methods to those introduced in Volume 2, Section 3.1), to the development of 
dedicated solvers which solve both the Navier-Stokes equations and neutron transport in one code. 


The methods used in the solution of the neutron transport equation are described below, ranked 
in an approximate order of their accuracy from higher fidelity steady solution methods to lower 
fidelity methods that can model the transient behaviour of reactors. There is a clear division in 
accuracy between high-fidelity probabilistic (random selection of events of known probability) and 
deterministic (no randomness) methods. 


Probabilistic Methods 


The tools that use the probabilistic approach are generally solved using Monte Carlo methods to 
numerically solve the neutron transport equation. The Monte Carlo approach involves the random 
selection of a range of probable events defined by the neutron cross sections, the mass of fuel, 
the presence of absorbers, moderators and reflectors and the temperature of each material in the 
reactor. These methods generally apply the continuous neutron cross sections. They can provide 
highly accurate solutions in ‘criticality searches’ when the behaviour of a large number (or popula- 
tion) of neutrons is solved to find the ke¢ of a particular reactor configuration. Accurate solutions 
are expected to have a 30 (where a is the standard deviation of k.~ obtained by the Monte Carlo 
solver) or 99.7 % confidence interval of around 10 pcm. The unit pcm is percent mille or 10—°, and 
is frequently encountered in reactor physics. 


The Monte Carlo solvers enable the user to approximate and homogenise the macroscopic cross 
sections into a number of discrete energy groups that can best approximate the neutronic be- 
haviour of the neutron spectrum of the reactor modelled. For example as few as two energy groups 
can be applied to some neutronics models that can be used to characterise the behaviour of an 
LWR, because the fast neutrons are quickly moderated to the thermal spectrum (CEA, 2015). In 
the case of reactor burn-up or for the modelling of fast spectrum reactors, more energy groups 
(hundreds) are required in order to approximate the capture of neutrons in the resonance region. 
In fast reactors, neutrons are not moderated by the coolant so they take much longer to reach the 
thermal or lower energy end of the neutron spectrum, hence more energy groups are required. 
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The results obtained by the Monte Carlo solvers can also be used as reference or benchmark 
solutions for lower fidelity deterministic methods. Common solvers include MCNP, MONK, Serpent, 
OpenMG, Shift and the KENO and TRITON packages from the SCALE" code system. 


Deterministic Methods 


Deterministic transport solvers involve the solution of the neutron transport equation discretised in 
the spatial domain, and may also involve the discretisation of the angular domain depending on 
the complexity of the model (increasing its accuracy and cost). The deterministic approach can be 
further sub-divided into: 


* Lattice-physics codes, which can condense and homogenise neutron cross sections. 


« Angularly discretised methods, diffusion approaches and point kinetics codes that need to 
use homogenised neutron cross sections as inputs. 


Lattice Codes 


Lattice physics codes are usually two-dimensional finely discretised (meshed) representations of 
sections of a reactor, resolving the contribution from individual components and materials of a fuel 
assembly or fuel pins. The solution methods either use the method of characteristics or collision 
probabilities to discretise the neutron transport equation. 


The Method of Characteristics (MOC) solves the ‘characteristic form’ of the neutron transport equa- 
tion (derived and expressed differently, but equivalent to the equation in Section 2.6.2, as described 
by Demaziére, 2020) by following straight neutron paths across a domain. Solvers include MPACT, 
OpenMOC, APOLLO2, CASMO5, WIMS/CACTUSOT, SCALE/NEWT and DRAGON. 


The collision probabilities method involves the spatial and angular discretisation of the ‘integral 
form’ of the neutron transport equation. The method uses ray-tracing to evaluate the collision prob- 
abilities to allow the discretised transport equation to be solved. Solvers include HELIOS2 and 
WIMS/PRIZE. 


Angular Discretisation Methods 


In tools that apply the spherical harmonics or discrete ordinates methods, the spatial variables are 
discretised on a mesh, while the angular variables are discretised by an approximation: 


¢ In the case of spherical harmonics or the Py method, angular variables are resolved as a 
series of spherical functions that describe the variation of the parameter on the surface of a 
sphere, where in this case WN is the order of the expansion. 


« In the case of the discrete ordinates (S,) method, angular variables are resolved at discrete 
intervals with a discrete set of vectors for the direction of the neutron flux. 


Volume 6 


Solvers which can apply the Py or Sy methods include FETCH2, SCALE/DENOVO, WIMS/TWOTRAN, 


Griffin and GeN-Foam (as the simplified SP3 method). 


12 www.ornl.gov/scale/overview 


58 of 89 


3.4.2.3 


3.4.2.4 


3.4.3 


Methodologies 


Diffusion Approximation 


If scattering and sources of neutrons can be approximated as isotropic (no angular effects), then 
the simplest tool used for modelling the spatial variation of neutron flux is the diffusion approxima- 
tion (using Fick’s law). The neutron flux is modelled using a continuum transport equation that is 
similar to those used for fluid mechanics. 


Nodal diffusion neutron transport is a finite difference method that has been commonly imple- 
mented in solvers including DIF3D, DYN3D, Griffin'?, SIMULATE5, PARCS, WIMS/SNAP. The 
nodes represent the fuel assemblies of a solid fuelled reactor and therefore, they are typically 
arranged in the same shape as the cross-section of the assemblies forming the reactor core, with 
either hexagonal prism or hexahedral nodes. Each fuel assembly is modelled with a number of 
nodes axially, to be able to determine the power distribution in the vertical direction. Finite element 
and finite volume versions of the diffusion transport equation have also been implemented. 


Point Kinetics 


The simplest approach is the application of point kinetics, where the reactor is reduced to a point or 
a series of points, and an ordinary differential equation approximates the change in the averaged 
neutron density in the reactor with time. This assumes that there is a constant shape of neutron 
flux distribution. 


Point kinetics solvers are typically implemented in or coupled to all system codes (see Section 3.5). 
Examples of point kinetics solvers applied to molten salt reactors include RELAP, SPECTRA, 
TRANSFORM, MELCOR, LiCore. 


Transient and Coupled Analysis 


Neutronics and thermal hydraulics calculations are coupled routinely (to the extent that it is com- 
mon to see the code pairs stated as a single name, such as TRACE/PARCS or ATHLET/DYN3D), 
and the fidelity used for each tool will depend on the problem at hand and what the most suit- 
able approximations are. Both diffusion approximation and point kinetics methods are capable of 
modelling the transient behaviour of nuclear reactors over a period of minutes, to hours and days 
with reasonably low computational costs. They are also suited to design iterations where a large 
number of conditions need to be assessed repeatedly, and are suited to being coupled to thermal 
hydraulics codes. 


In relation to MSRs, some recent examples of applying higher fidelity methods than system codes 
and point kinetics include: 
* Coupling DYN3D to CFD (Cartland Glover et al., 2019). 


* The GeN-Foam solver implemented in OpenFOAM, which solves the neutron diffusion or 
SP3 equations within the finite volume CFD solver framework (Fiorina et a/., 2017). 


* The MPACT code coupled to CTF as part of VERA (Volume 1, Section 4.4.2) to model the 
MSRE (Graham et al., 2021). 


13 Tools often contain multiple solution methods within a common framework for nuclear data processing, as an example 
Griffin also contains Sj solutions capabilities. WIMS and SCALE are similar — their different modules for each solution 
method are noted. 
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* The Pronghorn porous medium CFD code was coupled to the Griffin neutronics code to 
assess a gas cooled pebble bed benchmark (Novak et al., 2021) and an MSFR benchmark 
(Abou-Jaoude et al., 2021). 


* OpenMC coupled to Nek5000 (solving LES CFD) was recently demonstrated (Merzari et al., 
2021), representing the resolution of the smallest length scales in a large pebble bed geom- 
etry. 


System Analysis 


Molten salt reactors of any type need system analysis for the same reasons as other reactors, 
such as for the prediction of long duration transients including the behaviour of active control and 
protection systems or passive cooling. Krepel et a/. (2014) noted three requirements for system 
codes that are peculiar to modelling MSRs (although some designs do not require all of them): 


¢ The transport of DNPs. 
* Heat generation in fluids from fission. 


¢ Graphite moderator structures that can be hotter than the fuel (heated by irradiation and 
cooled by the salt). 


The last of these means that there is likely to be a greater need for more detailed and three- 
dimensional predictions of solid temperatures, and so different and more sophisticated ‘heat struc- 
tures’ may be needed compared to those used in typical system code analysis (Volume 2, Sec- 
tion 3.2). This can also be extended to consider the need for sufficiently good representations of the 
core internal structures, vessel and external structures that will play a role in passive heat removal 
conditions. 


To enable the prediction of solid and fluid temperatures, as well as flow rates, a number of aspects 
of a system code’s capabilities that are particularly relevant to MSRs need to be considered: 


« The low Reynolds number flows that may be encountered will require laminar representa- 
tions of flow losses and heat transfer correlations, which depend more on specific channel 
geometries than turbulent flow correlations (which can usually be used with a universal hy- 
draulic diameter). In addition, flows with Reynolds numbers in the region of laminar-turbulent 
transition can be challenging to model, and can lead to oscillatory numerical behaviour. 


¢ Modelling of natural circulation flows, which requires particular attention to be paid to hy- 
drostatic components of pressure variation. This may require different numerical solution 
approaches, particularly with laminar flows, because the driving forces and pressure drops 
can both be very small. 


¢ In addition, the modelling of flow conditions that are affected by buoyancy (such as heat 
transfer and influencing laminar-turbulent transition) requires convective heat transfer mod- 
elling able to identify the appropriate flow conditions and apply correlations that can represent 
mixed convection (Volume 3, Section 2.2.1). 


¢ Convection correlations that include entrance or development effects and account for tem- 
perature varying viscosity (which will depend on the local heat transfer) may be required to 
be sufficiently accurate. 


60 of 89 


Methodologies 


* The modelling of radiative heat transfer between surfaces (such as in an external reactor 
vessel cooling configuration) and between surfaces and adjacent semi-transparent fluids. 


+ If species transport needs to be modelled (for DNPs, for example), then the action of physical 
turbulent mixing may need to be included and the presence of numerical diffusion/dispersion 
(caused by low-order numerics and coarse nodalisation) needs to be controlled. 


The availability and maturity of these functions and attributes depends on the tool and its pre- 
vious applications. For example, system codes developed for two-phase accident predictions for 
LWRs will not necessarily have focussed on implementing and validating modelling of these phe- 
nomena. The available and widely applied system code tools introduced in Volume 1, Section 4.4 
are therefore less mature for their applicability to MSRs. Molten salt reactor designs are, however, 
amenable to the same modelling approaches as other technologies, with some modification, as 
demonstrated by a selection of recent examples, where either existing tools have been adapted, 
or new tools written, for MSR applications: 


TRACE has been used to model the MSRE (Xun, 2016), and also the AHTR (Wang et al., 2015). 
The latter required salt models to be implemented in TRACE (Richard et al., 2014). More 
generally the suitability of the US NRC’s tools for modelling MSRs (and other non-LWR tech- 
nologies) is reviewed in US NRC (2020a). 


RELAP was extended by Shi et a/. (2016) to add DNP transport to the point kinetics model, and 
a model of the MSRE was built, and validation comparisons made. A model of the MSBR 
design was also built and used to demonstrate a number of flow and reactivity transients. 


MELCOR was conventionally an LWR severe accident code, but is being developed so that it can 
represent non-LWR reactors. This includes modelling reactor normal operating conditions 
to use as initial conditions and inputs to the modelling of MSR transients and accidents. 
This capability has been demonstrated using MSRE flow-reactivity transient experiment data 
(Beeny et a/., 2022). It forms part of the US NRC’s intended route for the evaluation of source 
terms and accident consequences (US NRC, 2020b). 


TRANSFORM is an open source modelling tool developed by ORNL intended to allow energy 
system models to be rapidly developed using a library of connected components, and its 
models are potentially less geometrically detailed than other system codes (more akin to 
process simulation models). A model of the MSDR has been demonstrated using it (Green- 
wood et al., 2020). 


SAM is asystem code designed specifically for single-phase advanced reactors, and is being used 
by Kairos Power as their system analysis code (Blandford et al., 2020) for their pebble bed 
core FHR. An MSRE model has also recently been created and compared to flow-reactivity 
transients (pump coast-down and startup, as well as natural circulation, Hua ef a/., 2022). 
Some of the issues related to numerical diffusion and low order numerics described above 
are discussed in the context of SAM by Hu (2017). 


SPECTRA is an established system code that has been applied to a wide range of reactor types, 
and is under development to add additional functionality needed for MSRs. As well as DNP 
transport, modelling of the behaviour of fission products including dissolved gaseous species 
and noble metals is included and compared to MSRE data (Roelofs and Stempniewicz, 
2021). The same MSRE flow-reactivity transients as were modelling in SAM have been 
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demonstrated and the two codes compared as a recent benchmarking exercise (Hua et al., 
2022). 


Flownex is acommercial system code applied in a number of industries, with functionality suitable 
for MSRs. It has been applied to the design of the IMSR"4 and to FHR flows (Andreades et al., 
2016). It has also been used in comparisons for gas cooled pebble beds (Novak et a/., 2021). 


TREND is anew code written to support TMSR development, and demonstrate comparisons with 
the MSRE and also the CIET test facility (Yu et al., 2021). 


LiCore is in the early stages of development to support the assessment of the MSFR under the 
SAMOSAFER programme (Laureau et al., 2021). 


The reliance on the MSRE as a starting point for many modelling developments is apparent. This 
is partly because it is the only real reactor that can be used as a benchmark, but also because of 
the extensive design and test documentation available from that project. 


Thermal Hydraulic Experiments and Validation 


This section discusses flow and heat transfer experiments relevant to MSRs that have been (or are 
being) conducted, and what kind of data can be expected to be available from experimental meth- 
ods to help with model validation. Guidance on the different types of test (such as distinguishing 
Separate Effect Tests (SETs) from Integral Effect Tests (IETs)), to validation and to scaling is given 
in Volume 1, Sections 4.6 and 4.3.2 and in Volume 4. The availability and accuracy of thermophys- 
ical properties is highly relevant, but these experiments are not discussed here — Section 2.3 and 
the references mentioned provide that context. 


Two types of experimental approach are discussed: those using simulant fluids, and those using 
molten salts. 


Simulant Fluid Experiments 


The MSRE programme used a full-size mockup to study flow distribution and heat transfer. This is 
described in detail by Kedl (1970), and was modelled in SAM (Leandro et al., 2019) as a stage in 
building an MSRE system code model. The mockup used water, or water with additives to increase 
viscosity to provide a better Reynolds number match. 


The ability to use low temperature water, or another easier to handle simulant fluid, is common 
practice for some types of experimental investigation, and compared to high temperature salts, al- 
lows for faster development, cost reductions and high-fidelity measurements to be made. It can also 
enable access for optical techniques to provide spatially resolved, ‘CFD grade’ validation data. An 
example of this is the use of a fluid with a refractive index matched to a transparent acrylic, so that 
Particle Image Velocimetry (PIV) can ‘see through’ the solid objects in pebble beds (Nguyen et al., 
2018) and rod bundles (Nguyen et al., 2020) to resolve the flow fields in the spaces between them. 
Using alternative fluids, which have different material properties to the prototypical application, re- 
quires a scaling analysis to be performed to match the non-dimensional number characteristics 


14 flownex.com/industries/nuclear-smr 
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between the test and real reactor as well as possible, and to quantify scaling distortions. Com- 
pared to the two-phase flow and heat transfer requirements of LWR tests, however, matching the 
behaviour for single-phase flows is more straightforward. 


Heat transfer experiments with alternative fluids have additional complications in attempting to rep- 
resent the dependence of Nu on surface heat flux via the temperature dependence of viscosity on 
temperature. Laboratory experiments may not be able to work with the high heat fluxes necessary 
to create sufficient temperature differences, or the extent of the temperature difference necessary 
for the property variation may distort the applicability of the simulant fluid. Mixtures of several fluids 
will be needed to match the required temperature-viscosity relationship, as well as other thermo- 
physical properties and the refractive index. 


It has been found that the heat transfer oil Dowtherm A can be used as a simulant for FLiBe, as 
discussed by Andreades et al. (2016) and Liu et a/. (2018). The properties of water, Dowtherm 
A and FLiBe are compared in Figure 3.2. Three example properties relevant to thermal hydraulic 
performance are shown, and each fluid is plotted over its own temperature range from T, to Ty, 
normalised to T’ = (T — T,)/(Tu — Tz). This shows that, over the range chosen, the Prandtl 
numbers of Dowtherm A and FLiBe are very well matched, and while that of water is lower, it does 
decline in a similar ratio with increasing temperature. The values of 6 and pc, for water and FLiBe 
are comparable. 


Dowtherm A was applied in the Compact Integral Effects Test (CIET) facility to study natural- 
circulation driven decay heat removal in FHRs (Zweibaum ef al., 2016), and scaling analysis can 
be used to design facilities with reduced size and power where the Re, Gr and Pr numbers can 
be matched simultaneously to a reactor (Zweibaum ef a/., 2020). Blandford et al. (2020) describes 
how both water and oil based systems are being developed to test the flow and heat transfer per- 
formance of the Kairos Power FHR design, prior to the construction of a molten salt based (non 
nuclear) larger ‘Engineering Test Unit’. This ETU has objectives more focussed around proving 
the ability to build equipment at prototypical scale and conditions, and the integrated process of 
operating the nuclear Structures, Systems and Components (SSCs) (including pumping salts and 
handling fuel for example). It is not intended for investigating and demonstrating more fundamental 
thermal hydraulic phenomena and validating models, which is more suited to scaled SET and IET 
test facilities. 


Molten Salt Experiments 


The MSRE project documentation is a source of a large range of information, including thermal 
hydraulic test data (as noted in Section 3.5), but it is not a source of sufficiently well resolved 
results to be used as CFD grade validation data. 


A number of contemporary molten salt flow and heat transfer facilities are being developed. While 
they are less able to provide detailed flow and precise heat transfer validation data'®, they have 
an important role to play in allowing organisations to ‘to acquire some technical experience on 
the handling and processing of molten salts.’ (Pioro, 2016, Chapter 7) — providing a technology 


15 Temperature measurements are complicated by heat losses from test fixtures at high temperature — the details of insulation 
and external convective and radiative heat transfer can be significant. The accuracy of material property data at elevated 
temperatures may also not be sufficient. 


63 of 89 


Methodologies 
30 g x10" gsi? 
Water 
Dowtherm A 8 


FLiBe 


25/7 


| 


Volumetric Heat Capacity, pc, (J/K m’) 
is 


ipo} 

oO 
: 

fo) 
: 
i) 


a 
T 
wo 
oa 


Prandtl Number, Pr 
a 


= 
oO 
T 
wo 
T 


25/7 


Volumetric Thermal Expansion, (3 (1/K 
he) a 


= 
T 
ie) 
T 


——a 


15 


Normalised temperature scale, T’ 
Water: 5 to 60°C + Dowtherm A: 45 to 105°C —-FLiBe: 520 to 700°C 


Figure 3.2: Comparison of some properties of water, Dowtherm A and FLiBe. 
Data sources: FLiBe: Romatoski and Hu (2017) Water: IAPWS (Wagner et a/., 2000); 
Dowtherm A Dow (2021). 


development route for valves, pumps, heaters, insulation and materials. A substantial part of the 
challenge of obtaining valuable data is using instrumentation appropriate for the high temperature 
salt environment; this is reviewed by Holcomb et al. (2018). 


There are several experimental facilities associated with the SAMOFAR and SAMOSAFER projects: 


FFFER: Forced Fluoride Flow for Experimental Research. A recirculating loop with gas injection 
and a bubble separator, testing the efficiency of gas removal from flowing salt. 

SWATH: Salt at WAII: Thermal ExcHanges. A channel geometry to study heat transfer. Laminar 
and turbulent flows can be created and solidification at walls can be studied. 

DYNASTY: DYnamics of NAtural circulation for molten SalT internallY heated. A high aspect ratio 
loop geometry with heated pipe walls to replicate some of the features of a salt with internal 
heat generation undergoing natural circulation. 


Rubiolo et al. (2017) describes FFFER and SWATH. Both SWATH and DYNASTY use water ex- 
periments for development and flow validation purposes as a precursor to salt operation. 


DYNASTY has similarities to other molten salt loops, such as those described by Srivastava et al. 
(2016), Srivastava et al. (2017), and the data available is generally limited to integral parameters 
such as flow rates and individual temperature measurements at a given location. While these are 
suitable for the purpose of assessing natural convection stability as well as overall circulation and 
heat transfer (useful for assessing and validating system code models, and also the situation for 
water based experiments with similar geometries, Wilson et al., 2022), it does not provide spatially 
resolved data. 


Salt loop geometries with spatially resolved temperature fields are beginning to be developed. 
Arora et al. (2021) describe how ‘optical fiber distributed temperature sensors’ (combined with 
point measurements by thermocouples) are able to measure the temperature distribution of FLINaK 
under forced convection conditions over the cross-section of a channel and along its length. Britsch 
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et al. (2019) describe a FLiBe loop operating under natural circulation, with laminar flow that uses 
similar distributed fiber measurements. 


Local salt velocity or flow rate measurements are difficult to make, but bulk flow rate measure- 
ments can be made using indirect methods (for example tracking the progress of a temperature 
perturbation) or external sensors (such as a transit time ultrasonic flow meter). 


Freezing within salt experiments, whether desired or inadvertent is a significant factor. Britsch et al. 
(2019) describes how frozen salt layers formed in cooler sections, affecting heat transfer. Some 
experiments are specifically aimed at the assessment and development of ‘freeze valves’ and the 
validation of modelling predictions of them (Jiang et a/., 2020) — Chisholm et a/. (2020) provides a 
review of the design and operational experience of these components from the MSRE. 
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Most of the modelling methods for thermal hydraulics (and coupled neutronics) phenomena that are 
necessary for the prediction of single-phase flow and heat transfer in MSRs are relatively mature 
and available for industrial application. The areas of future development that will broaden these 
methods to be able to perform a wider range of assessments with more fidelity and accuracy (or 
less uncertainty) are: 


* The prediction of material properties using atomistic simulations. 


* Broadening the range of coupled multiphysics that interact with thermal hydraulics to include 
chemistry and fuel performance. 


These topics are discussed below to conclude this volume. 


Property Derivation from Atomistic Models 


Where material properties are unavailable and are hard to measure, it may be possible to de- 
termine them using atomistic simulations — broadly termed Molecular Dynamics (MD). While this 
section is aimed at molten salt properties, itis also applicable to nuclear fuel performance and other 
disciplines involving material science and physical chemistry. The descriptions here are intended 
for thermal hydraulics engineers who may wish to access MD derived quantities. The intention is to 
introduce the relevant vocabulary, provide an understanding of the physics that can be modelled, 
and what models of various scales and costs can be expected to achieve. It is not expected that 
these simulations will be performed in an industrial setting in the near future — specialists normally 
found within academic or research organisations are required. However, the modelling discipline is 
active, broadly practised and mature, with long standing international collaborative organisations’, 
so the necessary skills should be obtainable. 


In most engineering thermal hydraulics simulations, the atomistic nature of matter does not need 
to be considered and the continuum mechanics approximation is valid. ‘Constitutive equations’ 
mathematically express the role of the material properties, which are themselves a mathematical 
expression of the net effect of the very small scale physical interatomic (or intermolecular) forces. 
It is unusual in fluid mechanics to need to consider a fluid in a non-continuum manner directly”, but 
of knowledge microscopic behaviour can be used to provide macroscopic information. 


MD solves Newton’s equations of motion on a per-atom basis, where the forces (or potential en- 
ergies) between them are described by ‘interatomic potentials’ that describe the attractive and 
repulsive interaction between atoms as a function of their distance. For molecules, bonds can be 


1 Such as CECAM, www.cecam.org. 
2 Nano-scale liquid geometries and rarefied gas dynamics, such as high speed re-entry flows, are two examples. 
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treated like flexible springs or rigid links that hold the atoms together, and individual atoms in one 
molecule ‘feel’ the effect of many other atoms in other adjacent molecules. Net charge, polarisation 
and electric fields, and the associated longer range forces between atoms produced can also be 
included. 


In ‘classical’ MD simulations, the functional forms and coefficients of the potentials are empirical, 
with values for coefficients chosen to best replicate the properties of materials. Forces are generally 
assumed to apply only between pairs of atoms (so excluding simultaneous multi-atom interactions). 


The interatomic potentials can also be described using quantum mechanics?. One method to im- 
plement this is known as Density Functional Theory (DFT), which models the simultaneously in- 
teracting electronic structure of many atoms in their condensed state (liquid or solid). Ab initio* 
calculations merge DFT (or another quantum mechanical method) with molecular dynamics, which 
enables ‘on-the-fly’ calculation of the forces, and simulates the resulting motion of atoms (lftimie 
et al., 2005). 


Classical MD offers a significant (orders-of-magnitude) reduction in complexity and computational 
effort compared to ab initio calculations®, and is still able to make a substantial range of predic- 
tions, but is dependent on the applicability and accuracy of the potentials used. For simple fluids 
(liquid mono-atomics) and covalent molecules (from water to polymers and large biomolecules) the 
potentials are well described and tested, and are applicable across a large range of conditions (for 
example, they can model from solidification to high temperatures without modification or tuning). 
They can be derived or validated using many different types of data (including X-ray and neutron 
diffraction/scattering for example). However, where there is a substantial influence from electrons 
or ions in a material, such as in metals, or ionic liquids like salts, accurate and transferable classical 
potentials are harder to determine. 


MD is a computationally expensive discipline, usually requiring access to a High Performance 
Computing (HPC) facility. However, there are a large number of established software packages 
available, which are typically open source. 


« For classical MD, LAMMPS has a material properties focus, DL_POLY is a general purpose 
code, and GROMACS and NAMD are more focussed on biomolecules, but also able to be 
used for material property assessment. Classical MD can be used to determine a material's 
viscosity, heat capacity, thermal conductivity, density (and 6), diffusion coefficient, electrical 
conductivity and surface tension. Methods for measuring properties from a simulation are 
divided into 

— equilibrium methods, where the properties are determined from statistical fluctuations 
in a homogeneous sample of a material, and 


— non-equilibrium (NEMD) methods, which directly apply gradients of velocity or temper- 
ature, and measure the stress or heat flux response. 
Choosing which method to use is usually a case of that which will give the most converged 
(lowest variance) answer with the minimum computational cost, which depends on the num- 
From Schrédinger’s equation, and the Born-Oppenheimer approximation, which is the separation of electronic and nuclear 
motion, arising from the fact that nuclei are much more massive, so are considered stationary relative to the motion of the 
electrons. 


From latin, meaning literally ‘from the beginning’, but in scientific contexts, meaning ‘from first principles’. 
5 As an approximate analogy, ab initio MD is to classical MD as DNS (or LES) is to RANS CFD. 
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ber of atoms simulated (thousands to millions) and the duration modelled (picoseconds to 
nanoseconds). 


¢ For ab initio MD, commonly applied packages include ABINIT, CASTEP, CP2K, CPMD and 
VASP, and can also determine redox potentials and optical/IR spectral properties as well as 
being able to include chemical interactions. It is also possible to use ab initio simulations 
to derive coefficients for classical MD potentials, as discussed in the context of MSRs salts 
below. 


Examples of Using MD for Molten Salts 


There are a relatively large number of examples of the application of MD to the prediction of molten 
salt properties that can be found in the academic literature. These include storage/solar salts and 
MSR coolants and fuel salts. In common with the assessment of the outputs from any modelling 
field, caution is needed to distinguish the literature that is most rigorous and well developed (in 
the same way as for CFD research). In particular, the assessment given for whether a predic- 
tion requires a classical or ab initio approach, and for the former especially, where the input data 
(potentials) come from are key determinants of utility. 


For example, Madden and Salanne (2016) provide an introduction to the opportunities for using 
MD to determine the properties of fuel bearing salts, and also indicate some of the more subtle 
limitations of what can be expected, for example: 


* It is possible for different ionic states of the same cation® to be present (U* vs U**). These 
cations can be treated as a separate species, but the effect of charge transfer will need to be 
accommodated by part of the potential. 


« The presence of covalent bonds is neglected by assuming that salts are pure halide melts in 
classical models. The presence of other negatively charged ions, such as oxygen may inval- 
idate this assumption because some long-lived covalent species (such as UO7*) will form. 
This would require the use of ab initio MD to treat correctly, which may be computationally 
too expensive to be practical. 


Ab initio approaches have been recently used directly: 


¢ To assess the effect on redox potential of solute ions in MSR relevant salts (Nam and Morgan, 
2015). 


* To predict the effect of various metal impurities in fluoride salts on the IR absorption spectrum 


for thermal radiation (Chaleff et a/., 2018). 


However, for transport property derivation, they are computationally too expensive to apply widely. 
There is a significant body of ongoing research where classical molecular dynamics models are 
used to determine properties, applying a ‘polarizable ion model’, which are potentials obtained with 
the assistance of ab initio calculations. 


The process of determining the constants for the polarizable ion model potentials is described by 
Heaton et al. (2006) and Salanne et al. (2012), and demonstrated in a number of studies, such as 


§ Cations are positively charged, such as U**, anions are negatively charged, such as F. 
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those of Salanne et a/. (2009), Dewan et al. (2013), Gheribi et al. (2014) and Li et al. (2020). These 
papers report confident viscosity, specific heat capacity and density predictions, but only Gheribi 
et al. (2014) reports thermal conductivity. All references agree that conductivity is the hardest 
property to determine from an MD simulation. 


A recent combination of polarizable ion model MD and computational thermodynamics (Section 2.4.5) 
is demonstrated by Smith et a/. (2020) such that 


...a comprehensive structural-thermodynamic model of the LiF—BeF, [FLiBe] system 
is therefore reported that brings together the output of the experimental thermodynamic 
data and MD simulations. 


Again, however, the challenge of predicting thermal conductivity is highlighted. 


Coupled Multiphysics for Chemistry and Fuel Performance 


Multiphysics coupling of thermal hydraulics with neutronics is relatively mature, and important in 
MSRs. The neutronic and fuel performance of the reactor is highly influenced by the constituents 
of the salt, as are the redox conditions, and more coupling between chemical processes and other 
phenomena is expected to be of benefit to best estimate analysis. 


Chemical species transport models for MSRs are divided into two categories, depending on their 
purpose: 


Burnup/depletion: Aufiero eft a/. (2013), Betzler et al. (2017) and Hu et al. (2021) demonstrate 
the modelling of the evolution of fissile fuel cycles from startup to equilibrium including fast 
and thermal (in Betzler et a/., 2017) spectrum reactor models. They include online refuelling 
and reprocessing as well as the removal of fission products (for both nuclear reactivity and 
chemical reasons). These models are intended to account for the mass present in a system, 
and the effect that the presence of species (that transmute to other elements through nuclear 
reactions) have on criticality and source terms over time periods of years. 


Physico-chemical transport: The behaviour of DNPs and fission products such as '*°Xe (with a 
large thermal neutron absorption cross section) affect neutronic performance, as do helium 
bubbles injected to enhance the transport and removal of gaseous fission products. The re- 
dox potential (Section 2.4.1) of the salt influences the solubility of gaseous and noble metals 
fission products and the corrosion of vessel structural materials. Methods that can account 
for the complex and coupled mass transfer, mass accountancy, thermodynamic and chem- 
ical reaction effects are under development. For example Walker ef a/. (2021) describe a 
proposed test reactor experiment for gas and noble metal experiments in a molten salt. The 
test arrangement is modelled using CTF (which had additional functionality developed for 
this purpose) to predict surface deposition of fission products and their transport in bubbles 
and off-gases. Walker and Ji (2021) review MSRE data on fission product transport (and 
also provides a source of further relevant references) and applies the same CTF modelling 
to interpret the noble metal transport and deposition seen in those experiments. In physico- 
chemical models, accounting for the number of atoms of each element present is not suffi- 
cient for all purposes — how they form compounds (speciation) can have a significant effect, 
but modelling this is challenging in fuel salts because of the number of elements present, and 


69 of 89 


Future Developments 


the limited knowledge that exists for some of their physical and chemical reactivity properties 
(for lanthanide and actinide species in particular). 


Much of this modelling is ‘standalone’; it is not common for the more detailed chemical aspects 
to be modelled alongside burnup. To provide a ‘digital twin’ of a reactor through its life, future 
analyses that are more coupled or integrated are desirable, accounting for the evolution of species 
and reactor performance on long timescales, in the same framework as thermal hydraulics and 
neutronics. 


The transport of tritium is an example of this. There are existing standalone codes for tritium migra- 
tion (Kim et a/., 2010, Ohashi et al., 2008), developed with an emphasis on hydrogen production 
by Very High Temperature Reactors (VHTRs). The functionality of these would be useful to be 
included in a system level model that predicted other species transport, to make a unified tool. 


Two recent examples demonstrate the beginning of a unified approach, each starting from one of 
the two categories of model discussed above: 


« In most burnup models, the representation of thermal hydraulic conditions is very simple, and 
spatial resolution of the temperature field or DNP distribution is not included. Hu ef al. (2021) 
demonstrates a model of the effect on the rate of fuel addition and criticality as influenced by 
coupled, spatially resolved thermal hydraulics, DNP transport, neutronics and fission product 
removal. The MSFR model demonstrated tracks mass inventories over years, so the mod- 
elling and numerical approaches used are required to be suited for long transients. 


Graham ef al. (2021) describes how CTF, with modifications to add multispecies modelling, 
coupled to MPACT for neutronics, was applied as part of VERA to model DNPs, gaseous fis- 
sion products and helium bubbles. The Thermochimica tool (Section 2.4.5) was also coupled 
to provide thermodynamic and solubility predictions for the species included. VERA, when 
applied previously to LWRs, is also able to model coupled thermal hydraulics, neutronics and 
depletion for the long term evolution of fissile material (Kochunas et a/., 2017). 
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6 Nomenclature 


Latin Symbols 
A Area, m? 
At Atwood number (At = (p1 — p2)/(e1 + p2)) 
Bi Biot number (Bi = hL/kg) 
Cp, Cy Specific heat at constant pressure or volume, J kg~ K~* 
dorD_ Diameter (D, = 4A;;/Pcs for hydraulic diameter), m 
f Darcy-Weisbach friction factor 
Fo Fourier number (Fo = a,t/L?) 
Gr Grashof number (Gr = gL3Ap/v?p = gL?BAT/v?, using the Boussinesq approxi- 
mation Ap/p ~ —BAT, where AT is often taken as Ty, — Ts,00) 
g Acceleration due to gravity, ms~? 
h_ Specific enthalpy, J kg~+, Heat Transfer Coefficient (HTC), Wm~? K—! or height, m 
I Radiative intensity, W m~? sr~! or Wm~? sr~+ wm? for a spectral density, where sr 
(steradian) is solid angle 
J Radiosity, W m~? 
k Thermal conductivity, Wm-! K-? 
L Length or wall thickness, m 
M_ Molar mass of a species, kg kmol~? 
Ma _ Mach number (Ma = U/a, where a is the speed of sound) 
n_ Refractive index 
Nu Nusselt Number (Nu = AL/k¢) 
p Perimeter, m 
P Pressure (P; = static pressure, Py = total pressure), Nm~? or Pa 
Pe  Péclet number (Pe = RePr) 
Pr Prandtl number (Pr = cpu/ ke) 
q_ Heat flux (rate of heat transfer per unit area, g = Q/A), Wm~? 
Q_ Rate of heat transfer, W 
r Radius, m 
R_ Gas constant (for a particular gas, R = R/M), Jkg~!K~! 
R Universal gas constant, 8314.5 J kmol~! K~1 
Rip Thermal resistance, K W~! 
Ra Rayleigh number (Ra = GrPr) 
Re Reynolds number (Re = pUL/,, or for an internal flow Re = WD,,/A.<h) 
Ri Richardson number (Ri = Gr/Re’) 
Sr Strouhal number (Sr = fL/U, where f is frequency) 


85 of 89 


Nomenclature S, 


«x j<< 


Volume 6 


Stanton number (St = Nu/RePr) 

Time, s 

Temperature (7, = static temperature, 77 = total temperature), K 
Wall friction velocity (u, = \/Tw/p), ms~+ 

Velocity, ms~? or thermal transmittance, Wm? K~! 
Specific volume, m* kg~? 

Volume, m? 

Mass flow rate, kg s~+ 
Wall distance, m 


Non-dimensional wall distance (y* = yu, /v) 


Greek Symbols 


Qn SS FE YM KF Hm R WR 


4 


oa 


Thermal diffusivity (a = k/pcp), m?s~? 


Volumetric thermal expansion coefficient (8 = —(1/p)(0e/OT)), K~? 
Ratio of specific heats (y = cp/c,) 

Emissivity or surface roughness height, m 

Absorption coefficient, m~? 

Wavelength, m 


Viscosity, kg m~1s~? 


Kinematic viscosity and momentum diffusivity (v = w/p), m?s~? 
Density, kg m~? 

Stefan Boltzmann constant, 5.67 x 10-8 W m-? K~* 

Shear stress, Nm~? 


Porosity or void fraction 


Subscripts and Modifications 


b 
cs 
f 


i 


wn 


8 


~ 


Bulk (mass-averaged) quantity 
Cross-sectional quantity 

Quantity relating to a fluid 

Quantity relating to a particular species 
Total (stagnation) quantity 

Turbulent quantity 

Static quantity or quantity relating to a solid 
Quantity relating to a wall or surface 
Quantity far from a wall or in free-stream 
Average quantity 

Molar quantity 

Varying quantity 
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7 Abbreviations 


AGR 
AHTR 
ANL 
ARE 
ARIS 
CALPHAD 
CFD 
CHT 
CIET 
CSP 
DFT 
DHR 
DNP 
DNS 
DO 
DRACS 
DYNASTY 
ENDF 
FFFER 
FHR 
GEM 
Gen IV 
HPC 
HTC 
HTGR 
IAEA 


MOC 


Advanced Gas-cooled Reactor 

Advanced High Temperature Reactor 
Argonne National Laboratory 

Aircraft Reactor Experiment 

Advanced Reactor Information System 
CALculation of PHAse Diagrams 
Computational Fluid Dynamics 

Conjugate Heat Transfer 

Compact Integral Effects Test 

Concentrating Solar Power 

Density Functional Theory 

Decay Heat Removal 

Delayed Neutron Precursor 

Direct Numerical Simulation 

Discrete Ordinates 

Direct Reactor Auxiliary Cooling System 
DYnamics of NAtural circulation for molten SalT internallY heated 
Evaluated Nuclear Data File 

Forced Fluoride Flow for Experimental Research 
Fluoride-Salt-Cooled High-Temperature Reactor 
Gibbs Energy Minimisation 

Generation IV 

High Performance Computing 

Heat Transfer Coefficient 

High Temperature Gas-cooled Reactor 
International Atomic Energy Agency 

Integral Effect Test 

Integral Molten Salt Reactor 

Infrared 

Joint Evaluated Fission and Fusion 

European Commission Joint Research Centre 
Large Eddy Simulation 

Liquid-Fluoride Thorium Reactor 

Liquid Metal-cooled Fast Reactor 

Light Water Reactor 

Minor Actinide 

Monte Carlo 

Molten Chloride Fast Reactor 

Molecular Dynamics 

Method of Characteristics 
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MOSART MOlten Salt Actinide Recycler and Transforming 

MSBR Molten Salt Breeder Reactor 

MSDR Molten Salt Demonstration Reactor 

MSFR Molten Salt Fast Reactor 

MSR Molten Salt Reactor 

MSRE Molten Salt Reactor Experiment 

NEA Nuclear Energy Agency 

OECD Organisation for Economic Co-operation and Development 

ORNL Oak Ridge National Laboratory 

OSTI Office of Scientific and Technical Information 

PIRT Phenomena Identification and Ranking Table 

PIV Particle Image Velocimetry 

PSI Paul Scherrer Institute 

RANS Reynolds-Averaged Navier-Stokes 

RCCS Reactor Cavity Cooling System 

RTE Radiative Transfer Equation 

RVACS Reactor Vessel Auxiliary Cooling System 

SA Sensitivity Analysis 

SAMOFAR Safety Assessment of the Molten Salt Fast Reactor 

SAMOSAFER _ Severe Accident Modeling and Safety Assessment for Fluid-fuel Energy Reac- 

tors 

SET Separate Effect Test 

SFR Sodium-cooled Fast Reactor 

SSC Structure, System and Component 

SSR-W Stable Salt Reactor — Wasteburner 

SWATH Salt at WAII: Thermal ExcHanges 

TMSR Thorium Molten Salt Reactor 

TRISO TRistructural-|SOtropic 

UQ Uncertainty Quantification 

US DOE United States Department of Energy 

US NRC United States Nuclear Regulatory Commission 

UV Ultraviolet 

VERA Virtual Environment for Reactor Applications 

VHTR Very High Temperature Reactor 
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